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Abstract 



Using a theoretical framework based on the next-to-leading order Q CD-improved effective 
Hamiltonian and a factorization Ansatz for the hadronic matrix elements of the four-quark 
operators, we reassess branching fractions in two-body non-leptonic decays B — > PP, PV, VV, 
involving the lowest lying light pseudoscalar (P) and vector {V) mesons in the standard model. 
We work out the parametric dependence of the decay rates making use of the currently avail- 
able information on the weak mixing matrix elements, form factors, decay constants and quark 
masses. Using the sensitivity of the decay rates on the effective number of colors, Nc, as a 
criterion of theoretical predictivity, we classify all the current-current (tree) and penguin tran- 
sitions in five different classes. The recently measured charmless two-body B —>■ PP decays 
{B^ K^ri',B^ — > K^r]',B^ K^7i^,B^ n^K'^ and charge conjugates) are dominated 
by the A^^^-stable QCD penguins (class-IV transitions) and their estimates are consistent with 
data. The measured charmless B PV {B^ — >■ ujK^ , B^ uh^) and B — >■ VV transition 
{B — i> (j)K*), on the other hand, belong to the penguin (class- V) and tree (class-Ill) transitions. 
The class-V penguin transitions are iVc-sensitive and/or involve large cancellations among com- 
peting amplitudes making their decay rates in general more difficult to predict. Some of these 
transitions may also receive significant contributions from annihilation and/or final state inter- 
actions. We propose a number of tests of the factorization framework in terms of the ratios 
of branching ratios for some selected B —>■ /11/12 decays involving light hadrons hi and /12, 
which depend only moderately on the form factors. We also propose a set of measurements 
to determine the effective coefficients of the current-current and QCD penguin operators. The 
potential impact of B ^ hih2 decays on the CKM phenomenology is emphasized by analyzing 
a number of decay rates in the factorization framework. 
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1 Introduction 



Recent measurements by the CLEO Collaboration of a number of decays of the type 
B — >• /11/12, where hi and /12 are light hadrons such as hih2 = 7nr,TiK,rj'K,LjK, have triggered 
considerable theoretical interest in understanding two-body non-leptonic B decays. These de- 
cays involve the so-called tree (current-current) b —>■ {u,c) and/or b —>■ s {ot b d) penguin 
amplitudes with, in general, both the QCD and electroweak penguins participating. The ap- 
propriate theoretical framework to study these decays is that of an effective theory based on 
the Wilson operator product expansion obtained by integrating out the heavy degrees of 
freedom, which in the standard model (SM) are the top quark and bosons. This effective 
theory allows to separate the short- and long-distance physics and one can implement the per- 
turbative QCD improvements systematically in this approach. Leading order corrections have 
been known for quite some time and in many cases this program has been completed up to 
and including the next-to-leading order corrections [§. Present QCD technology, however, does 
not allow to undertake a complete calculation of the exclusive non-leptonic decay rates from 
first principles, such as provided by the lattice-QCD approach, as this requires the knowledge 
of the hadronic matrix elements < hih2\'Heff\B >, where Tie// is an effective Hamiltonian con- 
sisting of the four-quark and magnetic moment operators. These are too complicated objects 
to be calculated with the current lattice-QCD methods. Hence, a certain amount of model 
building involving these hadronic matrix elements is at present unavoidable. 

The approach which has often been employed in non-leptonic heavy hadron decays is based 
on factorization With the factorization Ansatz, the matrix elements < hih2\'Heff\B > 

can be expressed as a product of two factors < hi\Ji\B >< /i2|<^2|0 >. The resulting matrix 
elements of the current operators Ji are theoretically more tractable and have been mostly 
calculated in well-defined theoretical frameworks, such as Lattice-QCD |p|-pJ]l, QCD sum rules 
|]T2|-|l5[| and potential models [3,||T6Hl8|; some are also available from data on semileptonic and 
leptonic decays [0. One can then make quantitative predictions in this framework taking into 
account the theoretical and experimental dispersion in the input parameters in the decay rates. 
Factorization holds in the limit that one ignores soft non-perturbative effects. The rationale of 



this lies in the phenomenon of color-transparency [20|, in which one expects intuitively that a 



pair of fast moving (energetic) quarks in a color-singlet state effectively decouples from long- 
wavelength gluons. In the decays B hih2, with typically E^-^^.^ ~ 0{mB/2), the energy of 
the quarks leaving the interaction is large and soft final state interactions should be small and 
hence factorization should be a good approximation. Final state interactions generated by hard 
gluon exchanges are, however, perturbatively calculable and can be included. Phenomenology 
of the factorization hypothesis in the decays B D^*^tt{p), B J/^K^*^ and related ones, 
involving the so-called current-current amplitudes, has been worked out and compared with 
the existing data with the tentative conclusion that data in these decays can be described in 
terms of two phenomenological parameters, ai and 02 |@], whose values seem to be universal 



The decays B ^ hih2 have been studied repeatedly in the factorization framework [22-26 



However, with the measurements of some of the B — > /11/12 decays []T],3, theoretical interest 
in this field has resurged. In particular, NLO-improved perturbative framework with updated 
phenomenological input has been used in a number of recent papers p7|^ll to study the CLEO 
data. We would like to take a closer look at the non-leptonic two-body decays B — > hih2, in 
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which QCD and/or electroweak penguins are expected to play a significant role. 

There are several theoretical issues involved in B ^ hih2 decays, which one does not 
encounter in the transitions B — > i?i/i2, where Hi is an open {D^*\ D^*^) or bound {J/'ip, rjc, Xc) 
charmed hadron, or in decays such as B ^ D^*'>D^*\ which are governed by the current-current 
(tree) amplitudes. In the case of the induced b ^ s and b —>■ d transitions, penguins play an 
important role. Of these penguin transitions, the ones involving the top-quark can be reliably 
calculated in perturbation theory as they represent genuine short- distance contributions. The 
rest of the penguins, which involve both the charm- and light-quarks, also have genuine short- 
distance contributions which can be calculated using perturbation theory. Their importance in 
the context of direct CP asymmetries has been emphasized repeatedly in the literature p^ , |33[| . 
However, in principle, such penguin amplitudes may also involve significant non-perturbative 
(long-distance) contributions. Arguments for an enhanced role of non-perturbative penguin 



effects have been advanced in the literature |24[. In simpler cases, such as the electromagnetic 



decays i? — >■ + 7 and B —> K* + 7, charm-penguins are likewise present and they introduce 
(and higher order) power corrections akin to the long-distance effects being discussed in 
non-leptonic decays. In these cases, one finds that the power corrections are negligible 

p5|-^. The same holds for the non- resonant B —>■ Xgi^i^ decays [0. The pattern of 
the 1/m^-corrections remains to be investigated systematically for non-leptonic b {s,d)qq 
decays. However, it is suggestive that the next-to-leading order QCD-improved framework 
based on factorization can explain most of the recent CLEO data without invoking a significant 
non-perturbative penguin contribution [ p^j28|] . With improved measurements, this aspect will 
surely be scrutinized much more quantitatively. 

A related issue is that of the current- current b —>■ ccs and b ccd transitions feeding into the 



b — > sqq and b dqq transitions, respectively, by (soft) final state interactions (FSI) [p8l-|42 



While in the oft-studied case of 5 — > Kn decays, these effects are not found to be overwhelming 
for decay rates, yet, in general, it is not difficult to imagine situations where FSI may yield the 
dominant contribution to a decay width. There are three ways in which the amplitude for a 
decay in the factorization approach can become small: (i) the effective coefficients of the various 
operators entering into specific decays are small refiecting either their intrinsic (perturbative) 
values, implying they are small for iVc = 3, or their iVc-sensitivity meaning that they are small 
for some phenomenologically relevant value of Nc, (ii) due to CKM-suppression, (iii) due to 
delicate cancellations among various competing Feynman diagrams, resulting into an amplitude 
which is effectively small. Using A*",,, the effective number of colors, as a variable parameter, it 
becomes immediately clear that some linear combinations of the effective coefficients entering 
in specific decays are particularly sensitive to Nc and they indeed become very small for certain 
values of = 1/Nc. This then implies that other contributions such as the ones coming from 
FSI and/or annihilation may become important. A good case to illustrate this is the decay 
B^ — »• K^K, whose decay rate may be enhanced by an order of magnitude due to FSI [0 
and/or annihilation |^ contributions. 

In this paper, we undertake a comprehensive study, within the factorization framework, of 
all the two-body decay modes of the type B PP, B PV and B — > VV where P{V) is a 
light pseudoscalar (vector) meson in the fiavor U (3) nonet. Concentrating on the lowest lying 0~ 
and 1~ mesons, there are some seventy-six (76) such decays (and an equal number involving the 
charge conjugate states). The branching ratios of these decays are found to vary over four orders 
of magnitude. We calculate their decay rates (branching ratios) and work out the most sensitive 
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parametric dependence of these quantities. In many cases the factorized amphtudes are small 
due to the reasons mentioned in the preceding paragraph. While this by itself does not imply 
an intrinsic inability to calculate, it becomes difficult to be confident if the rate is additionally 
unstable, requiring a good deal of theoretical fine tuning in the factorization approach. We list 
all such two-body decay modes here and caution about drawing too quantitative conclusions 
on their widths based on the factorized amplitudes alone. We think that the sensitivity of some 
of the effective coefficients on Nc and the fine tuning required in some amplitudes can be 
used as a criterion of predictivity of i? — >■ /11/12 decay rates in the factorization approach. The 
pattern of color-suppression in current-current amplitudes has been previously used to classify 
the Ai'c-sensitivity of these decays into three classes 0. We extend this to also include the 
penguin-dominated decays, which belong either to iVc-stable (class-IV) or A^^c-sensitive (class- 
V) decays. In addition, penguin-dominated decay amplitudes involving large cancellations are 
also included in class- V. All penguin-dominated B PP decays belong to class-IV. This class 
includes in particular the decays — >■ K~^n~, B^ K~^7]', B^ K^t]' and B~^ — > 7i~^K^, 
measured recently by the CLEO collaboration (here and in what follows, charge conjugate 
decays are implied). On the other hand, the recently measured B PV and B VV decay 
modes by CLEO [0] are in class- V {B^ ujK^ and B K*(f)) or tree-dominated class-Ill 
[B^ — »• o'TT"'"). Possibly some of these, and many more examples of class-V decays worked out 
by us here, indicate that the factorization-based approach is rather uncertain in these decays 
and one may have to develop more powerful methods to make theoretically stable predictions 
in this class. Factorization approach is expected to do a better job in accounting for class-IV 
decays - a claim which is pursued here and which is supported by present data. 

We propose tests of factorization in 5 — ^ hih2 decays through measuring a number of ratios 
of the branching ratios which depend only on the form factors but are otherwise insensitive to 
other parameters, such as the effective coefficients and hence Nc-, quark masses, QCD-scale 
parameter and CKM matrix elements. The residual model dependence of these ratios on the 
form factors is worked out in two representative cases: (i) the Bauer-Stech-Wirbel (BSW) model 
IP and (ii) a hybrid approach, based on Lattice-QCD/Light-cone QCD sum rules, specifically 
making use of the results obtained in the frameworks of lattice-QCD [0,0 and the Light- 



Cone QCD sum rules [jT2|,|T5|. The proposed ratios will test factorization and determine the 
form factors. 

A quantitative test of the factorization approach lies in a consistent determination of the 
effective coefficients Oj of this framework. The QCD perturbative contributions to Oj can be 
calculated in terms of the renormalized Wilson coefficients in the effective Hamiltonian gov- 
erning the decays B hih2. Then, there are non-perturbative contributions which have to 
be determined phenomenologically. Of these ai and 02 govern the current-current amplitudes 
and they should be determined in B —>■ hih2 decays without any prior prejudice. Four of 
the Oj's (03, ...,a6) govern the QCD-penguin amplitudes and four more (07, ...,aio) govern the 
electroweak-penguin amplitudes. We propose measurements of selected branching ratios (and 
their ratios) to determine the effective coefficients CLi, CL2, CI4 and Oe from the first six from data 
on B hih2 decays in the future. Since the Wilson coefficients of the electroweak penguin 
operators in the SM are rather small in magnitude (except for Cg), which in turn yield very 
small branching ratios for these decays, a determination of 07, ...,aio is a formidable proposi- 
tion. The coefficient Og can be determined and we propose several decays to measure this. We 
also list decay modes in which electroweak penguins (hence 07, Oio) do play a noticeable role, 
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and work out their corresponding branching ratios. 

Finally, we explore the potential impact of the B ^ hih2 decays on the phenomenology of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [^. Here, we discuss relations of the type put 



forward by Fleischer and Mannel |^ (see, also |^6|) involving the decay rates of —>■ K^tt~ 
and B^ K^n^ , which can be used to determine cos 7, where 7 is one of the angles of the 
CKM unitarity triangle, in terms of the ratio of the tree-to-penguin amplitudes z = T / P and 5, 
the strong phase shift difference involving these amplitudes. A bound on sin^ 7 can be obtained, 
assuming that there are just the tree and QCD-penguin amplitudes: 

R = r(E± ^ n^K^) = 1 - 2 ^ COS 7 COS (5 + > sin^ 7 . (1) 
From this, constraints on 7 of the form 

0° < 7 < 7o V 180° - 7o < 7 < 180° (2) 
follow, where 70 is the maximum value of 7, which are complementary to the ones from the 



CKM unitarity fits [0,^. There are similar relations involving the decays B PV and 
B — > VV , where P = 7t,K and V = p, K* . A determination of the angle 7, however, requires 
the knowledge of Zi and 5i in these processes. Also, the effect of the electroweak penguins has to 
be included. Having a definite model, whose consistency can be checked in a number of decays, 
one could determine (within a certain range) the values of Zi and 5i. Given data, this would 
allow us in turn to determine 7 in a number of two-body non-leptonic B decays. We draw 
inferences on the angle 7 based on existing data on i?, and in line with |2^, we show that the 



allowed values of 7 (or the CKM-Wolfenstein |^9[ parameters p and 7]) from this analysis are 



consistent with the ones following from the CKM unitarity fits. Similar analysis can be carried 
out for the decays B —>■ PP, PV,VV, where now P = 7r°,7r^ and V = p^iP^- Measurements 
of these decays and their ratios would allow to draw inferences on the angle a. We illustrate 
this in the context of our model. The other kind of relations discussed by us involve ratios of 
the decay rates dominated by the 6 — > s and b d penguins, respectively. As pointed out in 
ref. |50], these ratios can be used to determine the ratio of the CKM matrix elements |V^d/Vjs|. 



Since this CKM-ratio will, in principle, be measured also in B^ - B^ mixings and radiative 
and semileptonic rare B decays |4^j5l|] , one could check the consistency of such determinations 
to reach quantitative conclusions about the QCD dynamics at work in non-leptonic decays. 
However, it is conceivable that some of the non-leptonic decays may already provide interesting 
information on Vtd before the other mentioned processes are actually measured. While not 
competitive in terms of eventual theoretical precision, non-leptonic decays are nevertheless 
quite instructive in this respect for the current CKM phenomenology. 

This paper is organized as follows: In section 2, we discuss the effective Hamiltonian together 
with the quark level matrix elements and the numerical values of the Wilson coefficients C^^^ in 
the effective Hamiltonian approach. In section 3, we introduce the factorization Ansatz, define 
the relevant matrix elements and discuss their evaluation in the BSW model and in the hybrid 
lattice QCD/QCD sum rule approach. The matrix elements for the three classes B PP, 
B — i> PV and B VV, obtained in the factorization approach, are relegated to Appendix A, 
B and C, respectively. Section 4 contains a discussion of the various input parameters (CKM 
matrix elements, quark masses, hadronic form factors and mesonic constants). The numerical 
input we use in the estimates of branching ratios are collected in various tables. In section 5, 
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we tabulate the values of the phenomenological parameters aj for three values of the effective 
number of colors (A^^^ = 2, 3, oo) for the four cases of interest b—^s,b^s,b^d and b ^ d. 
This serves to show both the relative magnitude of the effective coefficients of the various 
operators in i? — > hih2 decays in the factorization approach and also the stability of these 
coefficients against Nc- The classification of the B —>■ hih2 decays is also discussed here. We 
also discuss the contribution of the annihilation amplitudes and list some decays of potential 
interest. Section 6 contains the numerical results for the branching ratios which we tabulate for 
three specific values of the effective number of colors Nc = 2, 3, oo. The parametric dependence 
on ^ = 1/Ncis shown for some representative cases in various figures and compared with data, 
whenever available. In section 7, we list a number of ratios of branching ratios to test the 
hypothesis of factorization and give their values for the two sets of form factors (in the BSW 
and the hybrid Lattice-QCD/QCD sum rule approaches). We also discuss the determination of 
the effective coefficients ai, ag here through a number of relations. We estimate these ratios 
and make comparisons with data, whenever available. Potential impact of the B —>■ hih2 decay 
rates on CKM phenomenology are also discussed here. Finally, we conclude in section 8 with 
a summary and outlook. 



2 The Effective Hamiltonian 



2.1 Short-distance QCD corrections 

We write the effective Hamiltonian H^jf for the AB = 1 transitions as 



Gp 
V2 



10 



\i=3 



(3) 

where q = d,s and Ci are the Wilson coefficients evaluated at the renormalization scale fi. We 
specify below the operators in T^e// for 6 s transitions (for b ^ d transitions, one has to 
make the replacement s ^ d): 



0]* 


= SaYLUa ■ Up-fi^Lb/s , 


01 


= SaYLUf^ ■ 


U(3l,iLba , 




= Sa^LCa ■ Cp-f^Lbp , 


0% 


= Sal^Lci3 ■ 


Cp'^^Lba , 


03 


= Sal^Lba ■ Eg' Qpl p^Lq'p , 

= Sai'^Lba ■ J2q' q'pif^Rq'f^ , 


04 


= SaYLbp ■ 


Eg' Qpli^Lq'^ , 

Eg' q'pif.Rq'c, , 


05 




= Sc^YLbp ■ 


07 


= |s„7^L6„ ■ Eg' (iq'q'p^f,Lq'i3 , 


Og 


= \saYLbp 


■ Eg' eq'^fili,R(io, 

■ Eg' eq'q'p-i^^Lq'c^ 


09 


Oio 


= |s„7^L6^ 


0, 


= ((7,/87r2) ^ bp G% . 









(4) 



Here a and [3 are the SU{?>) color indices and A^^, A = 1, ...,8 are the Gell-Mann matrices; L 
and R are the left- and right-handed projection operators with L{R) = 1 — 75 (1 + 75), and 
denotes the gluonic field strength tensor. The sum over q' runs over the quark fields that 
are active at the scale /i = 0(mf,), i.e., (g'e{M, s, c, 6}). The usual tree-level IV-exchange 
contribution in the effective theory corresponds to Oi (with Ci{M\y) = 1 + 0{as)) and O2 
emerges due to the QCD corrections. The operators O3, . . . , Og arise from the QCD-penguin 
diagrams which contribute in order through the initial values of the Wilson coefficients at 
/i ~ My/ and operator mixing due to the QCD corrections |E3|. Similarly, the operators 
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Oj, . . . , Olo arise from the electroweak-penguin diagrams. Note that we neglect the effects of 
the electromagnetic penguin operator which we did not list explicitly. The effect of the weak 
annihilation and exchange diagrams will be discussed later. 

The renormalization group evolution from /i k, My/ to jU ~ mj, has been evaluated in 
leading order in the electromagnetic coupling and in the NLL precision in the strong coupling 
as [0]. Working consistently to the NLL precision, the coefficients Ci,...,Cio are needed in 
NLL precision, while it is sufficient to use the LL value for Cg. These coefficients depend on 
the renormalization scheme used. To obtain numerical values for the Cj we must specify the 
input parameters. We fix ^^(M^) = 0.118, aew{.Mz) = 1/128 and /i = 2.5 GeV. Then in the 
naive dimensional regularization (NDR) scheme, we have: 



(5) 



Here, C^'^'^ = Cg + C5. From the electroweak coefficients C7, ...,Cio, only Cg has a sizable 
value compared to the coefficients of the QCD penguins; its major contribution arises from 
the Z penguin. Note that the scale (/i) and scheme-dependence of the Wilson coefficients 
will cancel against the corresponding dependences in the matrix elements of the operators in 
H-effi as shown explicitly in |Q. Since, the matrix elements given below are obtained in the 



Ci 


= 1.117, 


C2 


= -0.257 , 


C3 


= 0.017 , 


C4 


= -0.044 , 


C5 


= 0.011 , 


Ce 


= -0.056 , 


C7 


= -1 X 10~^ , 


Cs 


= 5 X 10-^ 


Cg 


= -0.010 , 


Cio 


= 0.002 , 




= -0.158 . 







NDR-scheme, we have listed the values of the Wilson coefficients Q also in this scheme. 
2.2 Quark-level Matrix Elements 

In the NLL precision, the matrix elements of Tie// are to be treated at the one- loop level. The 
one-loop matrix elements can be rewritten in terms of the tree-level matrix elements of the 
effective operators 

(Vg'|K//|6) = E Cf^if^Ksq'mibr'- (6) 

id 

In the NDR renormalization scheme and for SU{3)c, the effective coefficients multiplying 
the matrix elements < sq'q'\Oj^^\b ■y^^^'^ become ( Vy and 7^ are the transpose of the matrices 
given below) 

c/^ = '^^ + ^y'v + iv log — j Q + ■ ■ ■ > 

C<)'' = C2 + — [ry + -fv log Cj + 



2i 



47r V " ° /i y 



Cr = C,-t^{Q + C, + Cg) + ^ rj; + ^^log'-^\ C,+ 



3i 



C'J^ = C,+ '-^{C, + C, + Cg) + ^ U + ^^log"^] C,+ 



247r ' ' ^ ^' 47r V /i / 



4j 



7 



cf 

CefT 
9 

'-"10 



1 a. 



a. 



^ Clew ^ 

Svr 



(7) 



We have separated the contributions Ct, Cp, and Cg arising from the penguin-type diagrams 
of the current-current operators Oi,2, the penguin-type diagrams of the operators O^-Oq, and 
the tree-level diagram of the dipole operator Og, respectively. Note also that we follow the 
procedure of ref. ||2^ of including the tree-level diagrams b ^ sg ^ sq'q' associated with the 
operator Og into the contribution Cg appearing in the expressions for C^^^ . So, we have the 
hadronic matrix elements of four-quark operators only. The process-independent contributions 
from the vertex- type diagrams are contained in the matrices ry and 7y. Here 7y is that part of 
the anomalous matrix which is due to the vertex (and self-energy) corrections. This part can 
be easily extracted from 7*^°) in ref. 



Iv 



( -2 
6 







V 
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-2 










-2 
6 







6 
-2 
2 




\ 




-6 
-16 / 



(8) 



The matrix ry contains constant, i.e., momentum- independent, parts associated with the vertex 
diagrams. This matrix can be extracted from the matrix f defined in eqn. (2.12) (and given 
explicitly in eqn. (4.6)) by Buras et al. in ref. [Q: 



rv 
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7 
3 


-7 














\ 




-7 


7 
3 
























63 
27 


-7 


















-7 


7 
3 
























1 

3 

-3 


1 




V 














35 
3 


/ 



(9) 



Note that the dependence and the scheme dependence of the vertex correction diagrams are 
fully taken into account in eqn. (|^ by the terms involving the matrices 7y and ry, respectively. 
There are, however, still scheme-independent, process-specific terms omitted as indicated by 



the ellipses, and we refer to |27] for a discussion of these omitted terms in exclusive two-body 
B decays. 

by 



The quantities Cp, and Cg are given in the NDR scheme (after MS renormalization) 



Ct 



E 



+ 



2 ml 




(10) 



8 



Cp — C3 



+ -log^ + -log^-AFi — -AFii—. 



3 3 ° /i^ 3 * 



i=u,d,s,c,b - 

_2m, .f 



3 °^ /i2 ^ I m? 



(11) 
(12) 



with Cg^^ = Cg + C5. The function AFi{z) is defined as 



4 / c/xx(l — x) log [1 — — x) — ie] 




AF,{z) -- 

The corresponding electroweak coefficient Cg is given by 



Ce = --(3C2 + Ci) E 



K'bKl, /2 2, ml 



q'q 



,t:, ^i^v;; V3 3 ^, 



- + -ln^- AFi ^ 



(13) 



(14) 



Note that the quantities Ct and Ce depend on the CKM matrix elements. In addition, the 
coefficients C^^ depend on fc^, where k is the momentum transferred by the gluon, photon or Z 
to the quark-antiquark pair q'q' in 6 — qq'q' . In two-body decays any information on k"^ is lost in 
the factorization assumption. However, given a specific model for the momentum distributions 
of the quark-antiquark pair inside the hadron, the partonic distributions calculated here can be 



folded with this distribution, as, for example, has been done in [^]. Since, we are interested 
here in the decays B —>■ hih2, where hi, h2 are light mesons, it is not unreasonable to assume 
that this smearing will be very similar in all the decays being considered. In particular, (fc^) 
is expected to be comparable in these decays. However, the actual value of (fc^) is model 
dependent. From simple two-body kinematics or from the investigations in ref. one 
expects k"^ to be typically in the range 



K < I, 
— ~ k 



2 < 



ml 



(15) 



As we shall see later, branching ratios considered here are not sensitive to the value of k"^ if it 
is varied in a reasonable range. 



3 Factorization Ansatz for the hadronic matrix elements 
of the four-quark operators 

We have now to calculate the hadronic matrix elements of the type {hih2\0i\B), where Oi 
are the four-quark operators listed in the preceding section. These will be calculated in the 
factorization assumption, which in the present context has been explained in a number of papers 
(see, for example, ref. |2^). To recapitulate briefiy, the hadronic matrix elements involving four- 



quark operators are split into a product of two matrix elements of the generic type {hi\qh\B) 
and {h2\q'q'\Q) 1 where Fierz transformation is used so that the fiavor quantum numbers of the 
quark currents match those of the hadrons. Since fierzing yields operators which are in the color 
singlet-singlet and octet-octet forms, this procedure results in general in matrix elements which 
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have the right flavor quantum numbers but involve both the singlet-singlet and octet-octet 
operators. No direct experimental information is available on the latter. In the factorization 
approximation, one discards the color octet-octet piece and compensates this by introducing a 
phenomenological parameter which determines the strength of the singlet-singlet contribution, 
renormalizing it from its perturbative value. The hadronic matrix elements resulting from the 
factorization are calculated in a model or determined from data, if available. 

To set our notation and introduce some auxiliary quantities which we shall need for numeri- 
cal calculations, we illustrate the salient features of our framework below. When a pseudoscalar 
meson is a decay product, such as in the decay B — PP, there are additional contributions 
from the {V + A) penguin operator Oq and Og- After Fierz reordering and factorization they 
contribute terms which involve a matrix element of the quark-density operators between a pseu- 
doscalar meson and the vacuum. For Oq involving 6 — > s transition {in b ^ d transition s is 
replaced by d), for example, this is given by 

< PiP2\Oe\B >= -2^(< Pi\sRq\0 >< P2\qLb\B > +[Pi ^ P2]) • (16) 

Using the Dirac equation, the matrix elements entering here can be rewritten in terms of 
those involving the usual {V — A) currents, 

< PiP2\Oq\B >= i?[Pi, P2] < PiP2\04\B > +[Pi ^ P2] , (17) 

with 

P Pi, P2 = ——f^ . . 18 

[niq + ms){mb - niq) 

Here, and rrig are the current masses of the quarks in the mesons Pi and P2. The same 
relations work for Og. Finally, one arrives at the form 

<PlP2\neff\B> = < P,\f\0 >< P2\j^\B > 

+ z2<P2|/^|o><Pi|j;|p> , (19) 

where j/^ and are the corresponding (neutral or charged) V — A currents. The quantities Zi 
and Z2 involve the effective coefficients, CKM factors and Gp. The 0~ — > 0~ form factors are 
defined as follows: 



{Pi{pi)\qi,Lh\B{pB)) 



\Vb + V\)ix ^ii 



FAQ') + '^^^q,Fo{q% (20) 



where q = Pb — Pi- In order to cancel the poles at q^ = 0, we must impose the condition 

Pi(0) = Po(0). 
The pseudoscalar decay constants are defined as: 

{P{p)\qrLq'\Q)=tfpp''. (21) 
With this, we can write the required matrix element in its factorized form 

(PiP2|7ie//|i^) = ^^ (^Q + q) fpMl - ml)F,^-^^{ml) + (1 ^ 2). (22) 
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The dynamical details are coded in the quantities a,, which we define as 



a,^Ct" + ^Ctl( (^ = odd); a,^C:" + ^Ctl( = even), (23) 

where i runs from i = 1, 10. Thus, we see that there are ten such quantities. They depend 
on the SM-input parameters, including the CKM matrix elements. The non-factorizing contri- 
butions in the matrix elements {hih2\0i\B) are modeled by treating Nc as a phenomenological 
parameter. Note, that this is the only place where Nc is treated as a phenomenological param- 
eter. In particular, in the calculation of C^^^ , we have used the QCD value Nc = 3. Insisting 
that there are no non-factorization effects present amounts to setting A*",, = 3 in calculating 
ttj. This is also referred to as "naive factorization' and is known not to work in decays such 
as i? ^ (D, D*)(7f, p), J/iJjK^*^ III these decays only the coefficients ai and 02 are 



determined. Note that QCD does not demand the equality of ai and 02 from these decays and 
from the ones B —>■ hih2, though their values may come out to be close to each other. Hence, 
all the ten quantities a, should be treated as phenomenological parameters and fitted from data 
on B hih2 decays. 

Returning to the discussion of the hadronic matrix elements, we recall that when a vector 
meson is involved in a decay, such as in i? ^ PV and B VV decays, we need also the 
B V form factors, which are defined as follows: 



{niB + niv] 
+ M(PB+Py)M 2 %]\^ 




q I rriB + my 



2mv{e* ■ q) 



I- 



q,Ao{q'), (24) 



where q = Pb — Pv: and e* is the polarization vector of V. To cancel the poles at = 0, we 
must have 

2myAo(0) = (niB + my)Ai(O) - (m^ - my)A2(0). (25) 
The decay constants of the vector mesons are defined as follows: 

< V\qi^iq\0 >= fv myef,. (26) 

This completes the discussion of the factorization Ansatz. The various input parameters 
needed to do numerical calculations, including the form factors and meson decay constants, are 
discussed in the next section. 



4 Input parameters 

The matrix elements for the decay B ^ hih2 derived in the preceding section depend on the 
effective coefficients ai,...,aio, quark masses, various form factors, decay constants, the CKM 
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parameters, the renormalization scale /i and the QCD scale parameter A^. We have fixed 
using the central value of the present world average as{Mz) = 0.118 ± 0.003 [^. The scale fi 



is varied between n = mb and /i = mb/2, but due to the inclusion of the NLL expressions the 
dependence of the decay rates on is small and hence not pursued any further. To be specific, 
we use fi = 2.5 GeV in the following. The dependence on the rest of the parameters is more 
pronounced and we discuss them below giving the present status of these quantities. 

4.1 CKM Matrix Elements 



The CKM matrix will be expressed in terms of the Wolfenstein parameters A, A, p and r]. 

I X A\\p - IT,) \ 

VcKM ^ -A (27) 

\ A\^{1- p-iT]) -A\^ 1 / 

Since the first two are well-determined with A = 0.81 ±0.06, A = sin^c = 0.2205 ±0.0018 p|, 
we fix them to their central values. The other two are correlated and are found to lie (at 95% 
C.L.) in the range 0.25 < t] < 0.52 and -0.2 < p< 0.35 from the CKM unitarity fits [g^. We 
shall show the dependence of the decay rates on the parameters p and rj in the allowed domain. 
However, for illustrative purposes and if not stated otherwise, we shall use p = 0A2,ri = 0.34, 
which are the "best-fit" values from the CKM unitarity fits WWi 

4.2 Quark masses 

The quark masses enter our analysis in two different ways. First, they arise from the contribu- 
tions of the penguin loops in connection with the function AFi{k'^ / mf) . We treat the internal 
quark masses in these loops as constituent masses rather than current masses. For them we 
use the following (renormalization scale-independent) values: 

rub = 4.88 GeV, = 1.5 GeV, = 0.5 GeV, m„ = = 0.2 GeV . (28) 

Variation in a reasonable range of these parameters does not change the numerical results of 
the branching ratios in question, as also investigated in [^]. The value of rub is fixed to be the 



current quark mass value fnb{p = 2.5 GeV) = 4.88 GeV, given below. Second, the quark masses 
rUb, rris, rrid and m„ appear through the equations of motion when working out the (factorized) 
hadronic matrix elements. In this case, the quark masses should be interpreted as current 
masses. It is worthwhile to discuss the spread in the quark masses, as determined from various 
calculational techniques and experiment. The top quark mass is now known rather precisely 
WLl{'mt) = 168 ± 6 GeV. Typical uncertainty on the 6-quark mass 6{TrTb{p = 2.5 GeV)) = ±0.2 
GeV ||58|j59[] is also small. Likewise, the mass difference rub — rUc = (3.39 ± 0.06) GeV ||59 



IS 



well determined, which can be used to determine rric reasonably accurately for the calculations 
being done here. Hence, to the accuracy of the present framework, the uncertainties in the 
decay rates related to 6mt, 6mb and 6mc are small and ignored. 

Light quark mass ratios have been investigated in chiral perturbation theory and up- 



dated in [P, yielding: m^/m^ = 0.553 ± 0.043, m,/md = 18.9 ± 0.9, m,/m„ = 34.4 ± 3.7. 

•^The corresponding "best-fit" values obtained in Q p ~ 0.15 and r/ ~ 0.34 are very close to the ones being 
used here. 
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Table 1: Input values in numerical calculations. 



names 


Values 


^ S \ ' ' A J 


0.118 


U 


2.5 GeV 


A 


0.81 


A 


0.2205 


t{B+) 


1.62 ps 




1.56 ps 


mt{'mt) 


168 GeV 


mb(2.5Ge1/) 


4.88 GeV 




1.5 GeV 




122 MeV 




7.6 MeV 


m„(2.5Ger) 


4.2 MeV 



These ratios were converted into the quark masses by using the QCD sum rule estimates of 
the s-quark mass of the somewhat older vintage rn7(l GeV) = 175 ± 25 MeV, yielding 
m^il GeV) = 5.1 ± 0.9 MeV, nQil GeV) = 9.3 ± 1.4 MeV Improved estimates based on 
QCD sum rules have been reported during the last year, which include 0(a^)-perturbative im- 



provements , improved estimates of yielding A||^ ~ 380 MeV, and improvements in the 



estimates of the spectral functions lowering the s-quark mass. A contemporary repre- 

sentative values of the s— quark mass in the QCD sum rule approach is: rn7(l GeV) = 150 ±30 
MeV [H. 

The corresponding estimates in the quenched lattice-QCD approach have been recently 
reported in a number of papers ||66|-|68|. The lattice community likes to quote the light quark 
masses at the scale fi = 2 GeV, and in comparing them with the QCD-sum rule results, quoted 
above for 1 GeV, one should multiply the lattice numbers by a factor 1.3. Representative 
lattice-QCD values are: ml(2 GeV) = 100 ± 12 MeV H, rn7(2 GeV) = 130 ± 2 ± 18 MeV H, 



and m7(2 GeV) = 110 ± 20 ± 11 MeV [Q. The error due to unquenching is largely unknown 
and for a discussion of the given lattice-specific errors, we refer to the original literature. Taking 
the last of these values as fairly representative, one now has the central value m7(l GeV) ~ 140 
MeV with a typical error of ±25 MeV - in reasonably good agreement with the QCD sum rule 
estimates. Using rnb{fi = rrib) = 4.45 GeV from the central value in p3] and 



m7(l GeV) = 150 MeV , m{l GeV) = 9.3 MeV , m;;{l GeV) = 5.1 MeV , (29) 

from the discussion above, the corresponding values at the scale /i = 2.5 GeV used in our 
calculations are given in Table 1. 

Varying the light quark masses by ±20% yields variation of up to ±25% in some selected 
decay rates (such as in —>■ rj'K^ and — » rj'K^, as also noted in [^). While this 
dependence should be kept in mind in fitting the quantities from precise data, this is clearly 
not warranted by present data. Also, fitting the values of the quantities aj is not the aim of 
this paper. Hence, we shall fix all the current quark masses to their values in Table 1. 
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Table 2: Form factors at zero momentum transfer in the BSW model [Q. 



Decay 


1 — x'n 


V 


A, 


Ao 


An 


B 


0.33 










B ^ K 


0.38 










B ^7] 


0.145 










B^i 


0.135 










B^p 




0.33 


0.28 


0.28 


0.28 


B K* 




0.37 


0.33 


0.33 


0.32 


B ^LU 




0.33 


0.28 


0.28 


0.28 



Table 3: Values of pole masses in GeV. 



Current 


m(O-) 


m(l-) 


m(l+) 


m(0+) 


ub 


5.2789 


5.3248 


5.37 


5.73 


db 


5.2792 


5.3248 


5.37 


5.73 


sb 


5.3693 


5.41 


5.82 


5.89 



4.3 Form factors and hadronic coupling constants 

Finally, we discuss the numerical values of the form factors and coupling constants introduced 
in the previous section. Concerning the form factors, we shall use two different theoretical 
approaches. The first is based on the quark model due to Bauer, Stech and Wirbel |^, which 
has been found to be rather successful in accommodating data on a number of exclusive decays. 
In the BSW model, the meson-meson matrix elements of the currents are evaluated from the 
overlap integrals of the corresponding wave functions. The dependence of the form factors on 
the momentum transfer squared (which in B ^ hih2 decays equals the mass squared of 
the light meson) is modeled by a single-pole Ansatz. The values of the form factors in the 
transitions B ^ ir, B ^ K, B ^ r], B ^ t]' , B ^ p, B K* and B ^ lu, evaluated at 
= are given in Table 2. We assume ideal mixing for the {u, 0) complex. This amounts to 
using in the quark language (j) = ss and u = -^{uu + dd). Note, that to implement the rj-i]' 
mixing, we shall use the two-mixing-angle formalism proposed recently in [^,0 , in which one 
has: 

|?7) = cos6'8|r78) -sin6'o|r7o) , W) = sin6'8|?78) + cos6'o|r7o) • (30) 

Here, rjg and tjq are, respectively, the flavor SU (3)-octet and -singlet components. The relations 
for the pseudoscalar decay constants in this mixing formalism involving the axial- vector currents 
A^, and A° are: 

{o\Ai\v{p)) = tf^p,, miwip)) = zf^,p„ 

m'Mp)) = ^/.%' miiv'ip)) = . (31) 

The best-fit values of the [rj-r]') mixing parameters from yields: ^8 = —22.2°, 6q = —9.1°, 
/8 = 168 MeV, and /o = 157 MeV, which are used to calculate the decay rates in which rj 
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Table 4: Form factors at zero momentum transfer from Lattice QCD and Light-cone QCD sum 
rules. 



Decay 


Fi — Fq 


V 




A2 


^0 


B 




0.30 ±0.04 










B 


K\Vo\ 


0.35 ±0.05 










B 


r] (see text) 


0.13 ±0.02 










B 


rj' (see text) 


0.12 ±0.02 










B 






0.35 ±0.05 


0.27 ±0.04 


0.26 ±0.04 


0.30 ±0.05 


B - 






0.48 ±0.09 


0.35 ±0.07 


0.34 ±0.06 


0.39 ±0.10 


B 


([11] & SU(3)) 




0.35 ±0.05 


0.27 ±0.04 


0.26 ±0.04 


0.30 ±0.05 



and/or rj' are involved. In deriving the expressions for the decays involving t] and t]', we include 
the anomaly term in and the contributions of 6 ^ sgg — s(?7, r^') as calculated in 



Definitions of the various matrix elements can be seen in the appendix and we refer to 
for further discussions. The values of the input pole masses used in calculating the form factors 
are given in Table 3. However, in the decays B hih2, only small extrapolations from = 
are involved, hence the error due to the assumed Q^-dependence and/or the specific values for 
the pole masses is small. 

The second and more modern approach to calculating decay form factors is a hybrid ap- 
proach, in which often lattice-QCD estimates in the so-called heavy light mesons, calculated 
at high-Q^, are combined with the Q^-dependence following from the light-cone QCD sum rule 



analysis [|T2|,|T3|. We refer to |]10| for detailed discussions, compilation of the lattice-QCD anal- 
ysis and references to the literature, and quote here the results from the UKQCD analysis 
jni. For the S ^ TT form factor: Fi(0) = Fo(0) = 0.27 ± 0.11; for S ^ p form factors: 
V{0) = 0.35lg;g^, Ai{0) = 0.27t°om, ^(0) = 0.26lE]:|]i, and Ao{0) = 0.30t[;:|]i The results 
from an improved light-cone QCD sum rule calculation [|T3| for Fi{B —>■ n) = Fo{B ^ tt) and 



Fi{B — > K) = Fo(B K) are given in Table 4. The results for Fi[B rj) = Fq(B ^ rf) and 
Fi[B 7]') = Fq{B ^ 7]') are calculated from the B tt form factors from [|l5l taking into 
account additionally the {rj, rj') mixing, as discussed earlier and further detailed in Appendix 
A. The results for the B — > K* form factors have been obtained in the light-cone QCD sum 



rule in ref. |T^, which yield: 



^i(O) 



B- 



Ai(0)^-^* 

y(o)^-^ 



ViO) 



0.76 ± 0.05 , (32) 
0.73 ±0.05 , 



which, in turn, lead to the estimates Ai(0)^^^* = 0.35 ± 0.07 and V{0)^^^* = 0.48 ± 0.09. 
Assuming similar S'f/(3)-breaking in the remaining two form factors, and using the estimates 
for the corresponding form factors in i? — >■ p quoted above, one gets: ^2(0)^^'^* = 0.34 ± 0.06 
and Ao{0)^^^* = 0.39±0.10. The values from this hybrid approach are collected in Table 4. As 
for the form factors in the BSW model, we use a simple pole approximation for calculating the 
form factors at different from = 0. However, for the decays of interest, this extrapolation 
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Table 5: Values of decay constants in MeV. 





fK 


fs 


fo 


fc 
J rj 


fri' 


/p 


Jk* 




/</> 


133 


158 


168 


157 


-0.9 


-2.3 


210 


214 


195 


233 



is small and one does not expect any significant error from this source. For example, for the 
B ^ P form factors, using the parameterization of Fq^i{Q'^) given in eq. (12) of ref. [|l^, the 



resulting difference in the form factors is found to be less than 2%. 

The values for the pseudoscalar and vector decay constants are given in Table 5. The values 
for fuj, fx, fx* and coincide with the central values quoted in extracted from data on the 
electromagnetic decays of u and r decays, respectively |]19[. The decay constants f^,, f^,, 
and defined in the appendix A are obtained from the values for /q and /§, ^8 and 6*0 for the 
(?7, 1]') mixing, given earlier. The errors on the decay constants in Table 5 are small (typically 
(1 — 3)%), except on /^f-* and f^'^^ for which we use here the estimates from |2^ obtained using 
the QCD-anomaly method. These quantities have also been determined from the r]c — v' ~ V~ 
mixing formalism and radiative decays J/ip ^ {rjc, rj' , r])'y and the two-photon decay widths 
{Vc,v'jV) ~^ 77 in ref. [27| with results similar to the corresponding values obtained using the 
QCD-anomaly method PS|. For some recent determinations of these quantities, see also |72 



5 Effective coefficients and a classification oi B ^ /11/12 
decays 

5.1 Effective coefficients 

The effective coefficients Oj, which are specific to the factorization approach, are the quantities 
of principal phenomeno logical interest. Note that there are four types of transitions that one 
encounters in the current-current and penguin- induced decays B ^ hih2: b s [b s], 
and b —>■ d \b —>■ d]. Numerical values of [i = 1,...,10) for representative values of the 
phenomenological parameter Nc are displayed in Table 6 and Table 7 for the 6 — *• s [6 — *• s] and 

6 — > (i [6 — > d] cases, respectively. A number of remarks on the entries in these tables is helpful 
for a discussion of the branching ratios worked out later. 

• The determination of oi and 02 in the 6 — c current- current transitions has received a 
lot of attention. It remains an open and interesting question if oi and 02 in the b ^ u 
transitions are close to their b ^ c counterparts, which have the phenomenological values 
Oi ~ 1 and 02 — 0.2 [T^J21[]. These values correspond to the parameter ^ = l/N^ having a 
value around 0.4. The decays B tttt, B ^ pir and B ujtt are well suited to determine 
these coefficients. 

• The coefficients 03 and in the QCD-penguin sector are smaller compared to 04 and ag. 
In particular, the combination 03 + 05 has a perturbative value of 3 x 10~^, i.e., for A^^ = 3, 
in all four cases resulting from large cancellations between 03 and 05. This coefficient also 
shows extreme sensitivity to the parameter Nc, which in the present model is a measure 
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Table 6: Numerical values of effective coefficients for 6 ^ s [ & ^ s] at A^c = 2, 3, oo, where 
Nc = oo corresponds to C^-^^ . The penguin coefficients Co/-^ . .... C^^^ and C^^^ are calculated 
for the Wolfenstein parameters p — 0.12 and r] — 0.34. Note that the entries for a3,...,aio have 
to be multiplied with 10~^. 









= 2 








= 3 








— oo 


ai 




0.99 


[0.99] 






1.05 


[1.05] 






1.16 


[1.16] 






0.25 


[0.25] 






0.053 


[0.053 






-0.33 


[-0.33] 


as 


-37 


- 14i 


[-36- 


- Ui] 




48 


[48] 




218 


+ 29i 


[215 + 29i] 


a4 


-402 


- 72i 


[-395 


- 72i] 


-439 


- 77i 


[-431 


- 77i] 


-511 


-87i 


[-503 - 87i] 




-150 


- 14i 


[-149 


- U{\ 




-45 


[-45] 




165 


+ 29i 


[162 + 29^] 


oe 


-547 


-72i 


[-541 


-72i] 


-575 


-77i 


[-568 


-77i] 


-630 


-S7i 


[-622 - 87i] 


a-j 


1.3- 


- 1.3i 


[1.4 - 


1.3i] 


0.5 - 


- 1.3i 


[0.5 - 


1.3i] 


-1.2 - 


- 1.3i 


[-1.1 - 1.3i] 


as 


4.4- 


- 0.7?: 


[4.4- 


0.7i] 


4.6 - 


- 0.4i 


[4.6 - 


0.4i] 




5.0 


[5.0] 


ag 


-91 - 


- 1.32 


[-91 - 


- 


-94^ 


- 1.3i 


[-94- 


- 1.3i] 


-101 ■ 


- 1.3i 


[-101 - 1.3i] 


a 10 


-31 - 


-0.7i 


[-31 - 


- 0.7i\ 


-14- 


-OAi 


[-14- 


- 0.4i] 




20 


[20] 



of non-factorizing effects. Hence, for decays whose decay widths depend dominantly on 
these coefficients, the factorization framework is not reliable. The reason is simply that 
the neglected contributions, such as the weak annihilation diagrams and/or feed down 
from final state interactions to these channels, could easily overwhelm the perturbative 
factorizable contributions. 

• Concerning the effective coefficients of the electrowcak operators, we note that 07, ag 
and oio are numerically very small. This again reflects their perturbative magnitudes, 
i.e. the coefficients C^'^'^ , as can be seen in the columns for Nc — 3. Varying one 
sees no noticeable enhancement in these coefficients (except for oio but it remains phe- 
nomenologically small to have any measurable effect). Hence, electroweak penguins enter 
dominantly through the operator O9, barring rather drastic enhancements (of O(IOO)) in 
the matrix elements of the operators O7, Os and Oio, which we discount. No attempts will 
be made to determine these coefficients here. In fact, in the context of the SM one could 
as well work with a much reduced basis in the effective theory in which the coefficients 
a/, ag and a 10 are set to zero. 

• The dominant coefficients are then ai, 02 (current-current amplitudes), 04, Gq (QCD 
penguins) and ag (electroweak penguin) , which can be eventually determined from exper- 
iments and we discuss this programmatically later. Of these ai, a2 (and to a very high 
accuracy also ag) do not depend on the CKM matrix elements. The dependence of a4 and 
oe (likewise, the smaller parameters 03 and 05) on the CKM factors enters through the 
function Cf. The numbers given in the tables for are obtained for the CKM parameters 
having the values p = 0.12 and rj — 0.34. Note that 02 depends strongly on Nc- 

This sets the stage for discussing the various branching ratios numerically and comparison 
with the available data. 
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Table 7: Numerical values of effective coefficients for h ^ d [h ^ d\ at A^^ = 2, 3, oo, where 
Nc = oo corresponds to C^^^ . The penguin coefficients C^^^ , ...,Cj^^ and C^^^ are calculated 
for the Wolfenstein parameters p = 0.12 and t] = 0.34. Note that the entries for a3,...,aio have 
to be multiphed with 10~^. 







iVe = 3 




ai 




0.99 [0.99] 






1.05 [1.05] 


1.16 [1.16] 


(32 




0.25 [0.25] 






0.053 [0.053] 


-0.33 [-0.33] 


03 


-33 


- 7i [-42 - 


■23i] 




48 [48] 


208 + IM [226 + 47i] 


04 


-377- 


34i [-423 


- im\ 


-412 


- 36i [-461 - 124i] 


-481 - Mi [-536 - 140i] 


as 


-145 - 


- Ui [-154 


- 14i] 




-45 [-45] 


155 + 14i [173 + 47i] 


ae 


-523 - 


3Ai [-568 


- 116i] 


-548 


- mi [-597 - 124i] 


-600 - 4H [-655 - 140z] 


ay 


1.5- 


l.Oi [1.1 - 


1.8i] 


0.7 


- l.Oi [0.3- 1.8i] 


-1.0 - l.Oi [-1.4 - 1.8i] 


ag 


4.5- 


0.5i [4.3 - 


0.9i] 


4.7 


- 0.3i [4.5 - 0.6i] 


5.0 [5.0] 


ag 


-91 - 


l.Oi [-91 - 


- 1.8z] 


-94 


- l.Oi [-95 - 1.8i] 


-101 - l.Oz [-101 - 1.8i] 


aio 


-30- 


0.5i [-31 - 


- 0.9i] 


-14 


-0.3z [-14-0.6i] 


20 [20] 



Before discussing the numerical results and their detailed comparison with experiment and 
existing results in the literature, it is worthwhile to organize the decays B ^ hih2 in terms of 
their sensitivity on and anticipated contributions due to the annihilation diagrams in some 
of these decays. 

5.2 Classification of factorized amplitudes 

In the context of the tree (T) decays, a classification was introduced in [Q, which is used 
widely in the literature in the analysis of B decays involving charmed hadrons. These classes, 
concentrating now on the B ^ hih2 decays, are the following 

• Class-I decays, involving those decays in which only a charged meson can be generated 
directly from a singlet current, as in B^ vr^vr", and the relevant coefficient for these 
decays is ai. This coefficient is stable against variation of Nc (see Tables 6 and 7). There 
are just five class-I decays: B^ 7r^7r+, B^ p 7r+, B^ P^^r^, B^ P P^i and 
exceptionally also B^ p^K^. 

• Class-II decays, involving those transitions in which the meson generated directly from 
the current is a neutral meson, like B^ tt^tt", and the relevant coefficient for these 
decays is 02, which shows a strong iVc-dependence (see Tables 6 and 7). There are twelve 
such decays 5° —>■ h\h\., where h\ and h\ are mesons from the set ir^ ,ri,ri' , p^ and u. The 
decays B^ n^r]^'^ exceptionally do not belong to this class, as their decay amplitudes 
proportional to a2 almost cancel due to the destructive interference in two tree diagrams 
having to do with the configuration vr^ ~ mm — dd and r]^'^ ~ {uu + dd) + .... Note that as 
a2 has the smallest value at Nc = 3, all class-II decays have their lowest values at Nc = 3. 

• Class-Ill decays, involving the interference of class-I and Class-II decays, as in this case 
both a charged and a neutral meson is present both of which can be generated through the 
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currents involved in i?e//- An example of these decays is 7r"'"7r°, and the relevant 

coefficient is ai +ra2, where r is process-dependent (but calculable in terms of the ratios 
of the form factors and decay constants). For r < 1, the A^"c- dependence of the class-Ill 
amplitudes is below ±20% w.r.t. the perturbative value. As we shall see, the quantity r 
may considerably enhance the A^c-dependence if r is well in excess of 1. This, in particular, 
is the case in 5+ — > p^Tr"*" and S+ — > lott'^ decays, where r ~ 2; hence these Class-Ill 
decays show marked Ai'c-dependence. However, one should note that the decay rates for 
this class do not have their minima at Nc = 3, but rather at Nc = oo, reflecting the 
behavior of ai -|- 02. There are eleven such decays involving — >• (tt"*", p"*") (tt^, 77, 77', a;) 
and exceptionally also the decay S+ — > K*~^r]', in which case the penguin amplitudes 
interfere destructively. Its decay rate is, however, rather stable w.r.t. the variation in Nc 
but small due to the CKM suppression. 

However, when QCD (P) and electroweak penguins (-Pgw) ^^^^ present, as is the case in 
the decays B — > hih2 being considered, in general, the above classification has to be extended. 
In this case, the generic decay amplitude depends on T + P + PgW- amplitude is still 

dominated by the tree amplitude, the BSW-classification given above can be applied as before. 
For those decays which are dominated by penguin amplitudes, i.e., T + P + Pg-\Y — P + -PgW' 
the above classification used for the tree amplitude is no longer applicable. 

For the penguin-dominated decays, we introduce two new classes: 

• Class-IV decays, consisting of decays whose amplitudes involve one (or more) of the 
dominant penguin coefficients 04, Qq and ag, with constructive interference among them. 
They are stable against variation in Nc (see tables 6 and 7) and have the generic form: 

M{B^ hfh^) ~ aioi + ^ a^ai + (33) 

1=4,6,9 

M(S°^/l?/l°) ~ 0:202+ 

j=4,6,9 

M{B^ ^ hfh2) ~ ai{ai+ra2)+ ^ ajOj + 

1=4,6,9 

with the second (P + Pg-yy) term dominant in each of the three amplitudes. The ellipses 
indicate possible contributions from the coefficients as, a^, 07, as and aio which can be 
neglected for this class of decays. The coefficients aj are process-dependent and contain 
the CKM matrix elements, form factors etc. The decays where cci and 0:2 are zero are pure 
penguin processes and are obviously included here. The tree-dominated decays, discussed 
earlier, also have a generic amplitude of the type shown above. However, in this case the 
penguin-related coefficients aj are numerically small due to the CKM factors (specifically 
due to Vtd < Vts). 

Examples of Class- IV decays are quite abundant. In our classification, all twelve B PP 
decays dominated by penguin amplitudes are class-lV decays. They include decays such as 
B^ — > ir+7r°, S+ — > K^r]^'\ which involve ai -|- ra2 as the tree amplitude, and S° — > X^tt", 
and S° K'^ri^'\ which involve 02 from the tree amplitude. Finally, the pure-penguin decays, 
such as B~^ — > n'^K'^, B^ — >■ K^K^ and P° K^K^ naturally belong here. There are 
altogether twenty nine such decays. The decay P° K*^r]', in contrast to its P "'"-counterpart. 
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is not a class-IV decay due to the destructive interference in the QCD-penguin amphtude. The 
variation in the decay rates belonging to class-IV decays is less than it30% compared to their 
perturbative {N^ = 3) value. 

• Class- V decays, involve penguins with strong Ai'c-dependent coefficients as, 05, 07 and aio, 
interfering significantly with one of the dominant penguin coefficients 04, ag and ag (anal- 
ogous to the class-Ill decays oi + ra2 dominated by tree amplitudes). Then, there are 
decays in which the dominant penguin coefficients (04, Gq, ag) interfere destructively. Their 
amplitudes can be written much like the ones in eq. (pSD, except that the sum in the sec- 
ond term now goes over all eight penguin coefficients. Since these amplitudes involve 
large and delicate cancellations, they are generally not stable against Nc- 

Examples of this class are present in B PV and B VV decays, such as B^ — > tt^4>, 
B° -> vrV, 5° r/W0, 5° -> uj(f), B^ p^(j), B° p°(j), etc. In all these cases, the 
amplitudes are proportional to the linear combination [03 + 05 — 1/2(07 + 09)] (see Appendix 
B and C). Examples of those where the amplitudes proportional to the dominant penguin 
coefficients interfere destructively are: B^ — B^ K^cj) etc. The above five classes 
exhaust all cases, though clearly there are some amplitudes where comparable T and penguin 
(P + P'Q^^) contributions are present. They can be assigned to one of the classes depending 
on their tree and/or penguin coefficients, the criterion being the A^^c-dependence of the decay 
rates. 

Summarizing the classification, Class-I and Class-IV decays are relatively large, unless sup- 
pressed by the CKM factors, and stable against variation of N^., which is a measure of non- 
factorizing effects in the present model. Class-Ill decays are mostly stable, except for the 
already mentioned B PV decays. Many Class-II and Class-V decays are rather unstable 
against variation of Nc either due the dependence on the iVc-sensitive coefficients or due to del- 
icate cancellations. Many decays in Class-II and Class-V may receive significant contribution 
from the annihilation diagrams which we discuss now. 



5.3 Contribution of annihilation amplitudes 

Annihilation (by which are meant here both W^^-exchange and ly^-annihilation) contributions 
are present in almost all decays of the type B hih2 being considered here. However, their 
contribution should be understood as power corrections in inverse powers of mt (equivalently 
in l/m^) in i? decays. In inclusive B decays, their contribution to the decay width relative to 
that of the parton model is determined by the factor 

mf \ rrih J 

where Jb — 200 MeV is the 5- meson decay constant. The near equality of the lifetimes of 
B^, B^{B^) and B^{B^) mesons shows that the above crude estimate is largely correct, and 
that annihilation contributions are sufficiently power-suppressed in 5-meson decays. For more 
sophisticated but in their spirit essentially similar calculations, see, for example, . 

However, in exclusive two-body 5-decays, the contribution to a particular channel depends 
on the CKM factors and the dynamical quantities Oj, and in some cases the non-annihilation 
contribution is enormously suppressed. In these channels, the annihilation diagrams, despite 
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being power suppressed in may yield the dominant contributions to the decay and must 

therefore be included in the rate estimates and CP-asymmctrics. Instead of working out the 
annihilation contribution in all the channels discussed here, which necessarily introduces un- 
known hadronic quantities, we do a classification of annihilation diagrams and list only those 
decays in which they are anticipated to be important. 

For the decays B — > hih2, we need to consider the following annihilation amplitudes: 

• Vr=^-Exchange: M{bd uu) =^ M{B^ {uq){qu)) oc 02 , 

• H^=^-Annihilation: M(hu du) =^ M{B^ {dq){qu)) oc aiA^ , 

• 1^=^- Annihilation: Mihu su) M{B'^ {sq){qu)) oc aiA^, 

where A = sin6'c- Here, qq is a light quark-antiquark pair. These amplitudes can be termed 
as the tree-annihilation contributions. In addition, there are also the penguin-annihilation 
contributions which are important for certain decays. For example, they feed dominantly to 

the decay B^ (fxf). 

There are yet more decays which can be reached via annihilation followed by rearrange- 
ment of the quark-antiquark pairs in the final state. Representative of these are the decays 
S± ^ (/.7r±, S± ^ and 50(5°) ^ ^tt^, ^^(SO) ^ (/.ryW, 50(5°) ^ (/.a; and 50(5°) ^ 
However, these rescattering effects (final state interactions) are expected to suffer from sup- 
pression due to the color-transparency argument used in defense of the factorization Ansatz. 
Since we have neglected these rescattering contributions in the factorization amplitudes worked 
out in this paper, it is only consistent that we also drop the annihilation contributions which 
feed into other channels through rescattering. 

We specify below those two-body B decays which are accessible directly in annihilation 
processes and hence may have significant annihilation contributions: 

• B^PP decays: 5° ^ Tr^r^W, 5° ^ r)r]' . 

• B ^ PV decays: B^ p^tt°, B^ p^ri^'\ 5° ^ a;7r°, B° ujr]^'\ B+ K*+K^, 
B+ K+(j), B° K*+K-, S° ^ K+K*-. 

• B ^ VV decays: B^ p^p^, B^ p^oj, B^ uu, B^ 00, B+ K*+K*^, 
B+ K*+<j), 50 ^ K*+K*-. 

Note, that in addition to the decay modes listed above, there are quite a few others in the Class- 
I, Class-Ill and Class-IV decays given in the tables, which also have annihilation contributions 
but in view of the large T and/or P + -PgW contributions in these decays, the annihilation 
contributions are not expected to alter the decay rates in these channels significantly and hence 
we have not listed them. 

The annihilation amplitude can be written as 

< hh^lHefflB Z < hh2\f\0 >< 0\j^\B > . (35) 

If hi and h2 are two pseudoscalars, the annihilation form factors are defined as 



111.^ III.. — _ //^^ /KVl2 



F^^m + ^^^^oP^Q^F.^'^'^m, (36) 
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where Q = Pi + P2- With this, we can write the required matrix element from the annihilation 
contribution (denoted here by a subscript) in its factorized form 



< PiP2\neff\B >a 



(37) 



where a^, i = 1,2. Note that the annihilation amplitude in the decay B —>■ P1P2 is proportional 
to the mass difference of the two mesons in the final state. Hence, in the present framework, 
there is no annihilation contribution to the decays such as B^ vr^vr^, B*^ — > K^K~ etc. 
Comparing this amplitude with the non-annihilation contributions given in eqn (^2]), one finds 
that the annihilation amphtude in 5 — > P1P2 decays is indeed suppressed by a hefty factor 



[m 



2 _ 



2)^-0 



(m 



B) 



[mi 



ml) 



(38) 



The annihilation form factors are difficult to relate directly to experimental measurements but 
they can be modeled. We expect F(f^^^(0) to have a similar magnitude as the the corresponding 



form factors F^^^^ (0), to which they are related by crossing, and which we have listed in Tables 



2 and 4. Based on this, the annihilation form factors appearing in eqs. (|37| ) and (^) are 
suppressed due to large momentum transfer at = m^, at which they have to be evaluated. 
The total suppression factor in B ^ PP decays is then (9(m^2/'^B)- However, the effective 
coefficients a^, i = 1,2 entering in the annihilation amplitude are much larger than aj, j = 
3, 10 governing the penguin-amplitudes. So, a part of the power suppression is offset by the 
favorable effective coefficients. 

In the decays B —>■ PV and B VV, we do not anticipate the annihilation suppression as 
severe as in the decay B —>■ PP. Concentrating on the decays B —>■ PV, the annihilation form 
factors are 



< PV\f\0 > 



—I 



mp + my 



(Pp - Pv),, - 
2my 



mp + mv)Ai{Q 



m. 



Q)Ao{Q') 



■a 



Q) 



A2{Q') 

mp + my 



(39) 



The annihilation matrix element in the factorization approximation can now be written as 
follows: 

< PV\neff\B >a= iV2GFV,bV;^,aaBmy{e* ■ pB)Ao(m|). (40) 

From this, it is easy to see that for this class of decays the suppression factor is only due to 
the large momentum transfer involved in the form factors Ao(m^). Hence, the annihilation 
diagrams can contribute more significantly in the decay amplitude. For some of the channels 
for which the non-annihilation contributions are highly suppressed, the annihilation diagram 
can be easily dominant. For example, the annihilation amplitude to the decay B^ — > K*~^K^ 
is 

ml). (41) 



< K*^K''\neff\B+ >a= tV2GFV:,Vu,ajBmK4e* ■pb)Ao( 
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If we take Ao{0) = 0.4, fs = 200 MeV, the annihilation branching ratio is of the order 10~^ 
which is an order of magnitude higher than the branching ratio calculated with the penguin 
contribution alone. Other channels where the annihilation channel may play a significant role 
have been listed above. 

For B VV decays, the conclusion is quite similar to the one for the B —>■ PV decays. 
However, as these decays involve yet more untested form factors, their numerical estimates 
require a model for these form factors. The suspected channels in B VV decays sensitive to 
annihilation contribution have been listed above. We conclude that the decays most sensitive to 
the annihilation channel are indeed the Class-II and Class- V decays, mostly involving B^{B^) 
decays. 

6 Branching Ratios and Comparison with Data 



Table 8: B ^ PP Branching Ratios (in units of 10"^) using the BSW [Lattice QCD/QCD sum 
rule] form factors, with k"^ = ml/2, p = 0.12, rj = 0.34, and Nc = 2,3,cxd in the factorization 
approach. The last column contains measured branching ratios and upper limits (90% C.L.) 
i- 



Channel 


Class 




= 2 


Nc 


= 3 




= 00 


Exp. 


5" 




I 


9.0 


[11] 


10.0 [12] 


12 


[15] 


< 15 


50 




II 


0.35 


[0.42] 


0.12 


[0.14] 


0.63 


[0.75] 


< 9.3 


B^ 


— i> 7]' 7]' 


II 


0.05 


[0.07] 


0.02 


[0.02] 


0.09 


[0.10] 


< 47 


B^ 


— > rjrj' 


II 


0.19 


[0.22] 


0.08 


[0.10] 


0.29 


[0.34] 


< 27 


B^ 


— > rjrj 


II 


0.17 


[0.20] 


0.10 


[0.11] 


0.24 


[0.29] 


< 18 


B+ 




III 


6.8 


[8.1] 


5.4 


[6.4] 


3.0 


[3.6] 


< 20 


B+ 


— s> Ti^rj' 


III 


2.7 


[3.2] 


2.1 


[2.5] 


1.1 


[1.4] 


< 31 


B+ 


71^1] 


III 


3.9 


[4.7] 


3.1 


[3.7] 


1.9 


[2.2] 


< 15 


B^ 




IV 


0.06 


[0.07] 


0.07 


[0.09] 


0.11 


[0.13] 


< 11 


B^ 


TT^T] 


IV 


0.20 


[0.24] 


0.23 


[0.27] 


0.30 


[0.36] 


< 8 


B+ 




IV 


9.4 


[11] 


10 


[12] 


12 


[15] 


< 16 


50 


K+TT~ 


IV 


14 


[16] 


15 


[18] 


18 


[21] 


151^ ± 1 


fiO 




IV 


5.0 


[5.9] 


5.7 


[6.8] 


7.4 


[8.9] 


< 41 


fi+ 


K+7]' 


IV 


21 


[25] 


25 


[29] 


35 


[41] 


65lii ± 9 


B^ 


K°7]' 


IV 


20 


[24] 


25 


[29] 


35 


[41] 


471^7 ± 9 


B+ 


K+7] 


IV 


2.0 


[2.3] 


2.4 


[2.7] 


3.4 


[3.9] 


< 14 


B^ 


K^T] 


IV 


1.7 


[1.9] 


2.0 


[2.2] 


2.6 


[3.0] 


< 33 


B+ 




IV 


14 


[17] 


16 


[20] 


22 


[26] 


23^iJ ± 4 


B+ 




IV 


0.82 


[0.95] 


0.96 


[1.1] 


1.3 


[1.5] 


< 21 


50 




IV 


0.79 


[0.92] 


0.92 [1.1] 


1.2 


[1.4] 


< 17 



The decay branching ratios are shown in Tables 8-11 for the decays B PP, B — *■ 
PV (involving h ^ d transitions), B PV (involving b —>■ s transitions) and B VV, 
respectively, for the two sets of form factors given in Tables 2 and 4. The numbers shown 
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for the hybrid Lattice-QCD/QCD sum rules correspond to using Ff^'^ = 0.36, F^^^ = 
0.41, Fi Q = 0.16 and F^ q = 0.145. The first two are shghtly above the range determined 



m 



T5| but within the (larger) range as determined from the lattice-QCD calculations [|Tl[]. This 
choice is dictated by data, as discussed in detail below. The fc^-dependence of the branching 
ratios in the range /c^ = ml/2 ± 2 GeV^ is small and hence the numbers in these tables are 
shown only for the case /c^ = ml/2. The CKM parameters are fixed at their "best-fit" values: 
p = 0.12,1] = 0.34. All other parameters have their central values, discussed in the preceding 
section. In these tables we give the averages of the branching fractions of 5° and B^, and of 
and B~ , respectively. Hence, when we refer to branching fractions in the following sections we 
always mean the averages over the B and anti--B decays. The CP-asymmetries are, however, 
in general quite sensitive to p3| , p2| . We shall discuss this point in a forthcoming paper on 
CP asymmetries [[f^ . 

A number of observations are in order: 

• There are so far five measured B ^ hih2 decay modes in well-identified final states: B'^ — >■ 
K+TT~, B+ K+T]', 5° K°7]', B+ TT+ir", and B+ ujK+, with their branching 
ratios (averaged over the charge conjugate modes) given in Tables 8 and 10. In addition, 
the decay modes B~^ — > Ti^h^iji^ = 7t~^,K~^) with a branching ratio B{B^ 7c^h~^) = 
{l.6toi ± 0.4) X 10-^ [g, the decay mode B+ Luh+{h+ = 7r+, ir+) with a branching 
ratio B{B+ -> ujh+) = (2.5t[!:?±0.3) x 10"^ and the decay modes B K*(j), averaged 
over B+ and S° decays with a branching ratio B{B K*(t)) = (l.ltoJ ± 0.2) x 10"^ @, 
have also been measured. 

• The branching ratios for B^ K^n^ and B^ n'^K'^ are in good agreement with the 
CLEO data. Moreover, being class-IV decays, they show only a small sensitivity on i^. 
The estimated branching ratios for B~^ — > tt+tt" and B~^ A'+vr^ are in agreement with 
the respective upper bounds. The latter being a class-IV decay is again stable w.r.t. the 
variation of N^, the former (a class-Ill decay) varies by approximately a factor 2.3 as Nc is 
varied. The branching ratio for the sum B^ ir^h^ is plotted as a function of ^ = l/N^ 
in Fig. |l| for the BSW model form factors (dashed-dotted curve) and two different sets, 
corresponding to the central values of the hybrid Lattice QCD/QCD-SR form factors 
(dashed curve) and for values which are closer to their theoretical range given in Table 4 
(dotted curve). We see that data for this mode is well explained. 

• We estimate the branching ratio for B^ —>■ tt+tt^ to be around 1 x 10~^ for the central 
values of the CKM parameters, which could go down to about 5 x 10~^ for Vub/Vct = 0.06. 
The present CLEO upper limit is in comfortable accord with our estimates but we expect 
that this decay mode should be measured soon. However, the decay B^ tt^tt^ is 
not expected to go above 10~^, which makes it at least a factor 10 below the present 
experimental sensitivity. 

• We show the dependence of the branching ratios on the input form factors and the param- 
eter ^ = 1/iVc for the decays B+ K+t]' and S° K° rj' in Figs. |^ and |^, respectively. 
As can be seen in these figures, data tends to prefer somewhat larger values for the form 
factors Fi o than the central values given by the Lattice-QCD/QCD sum rules in Table 4. 
However, the experimentally preferred values of the form factors all lie within the range 
allowed by the present theoretical estimates. Likewise, the branching ratio increases as 
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Table 9: -B ^ PV Branching Ratios (in units of 10~^) involving b —>■ d (or AS = 0) transitions 
using the BSW [Lattice QCD/QCD sum rule] form factors, with k"^ = ml/2, p = 0.12, r] = 0.34, 
and Nc = 2, 3, cx) in the factorization approach. The last column contains upper limits (90% 



CL.) from 
tables [|T9l. 



M. The upper limit on the branching ratio for p+vr^ is taken from the PDG 



Channel 



Class 



= 2 

q:q\ 

25] 

0.88] 

0.33] 

0.03] 

0.01] 

0.54] 

0.34] 

7.3] 

16] 

7.9] 

7.4] 

5.3] 

0.36] 

0.37] 

0.002] 

0.047] 

0.023] 

0.010] 

0.007] 

0.002] 



[28] 
[0.08] 
[0.10] 
[0.02] 
[0.001] 
[0.06] 
[0.02] 
[4.5] 
[15] 
[4.9] 
[6.5] 
[4.7] 
[0.44] 
[0.46] 
[0.0007] 
[0.005] 
[0.003] 
[0.001] 
[0.0008] 
[0.0006] 



Nr 



oo 
W0\ 
34] 
1.7] 
0.12] 
0.07] 
0.04] 
0.74] 
0.54] 
0.98] 
13] 
1.1] 
5.0] 
3.7] 
0.64] 
0.67] 
0.002] 
0.43] 
0.21] 
0.087] 
0.064] 
0.002] 



fio - 

B'- 

B^- 

B^- 

5° ■ 

fio - 

B^- 

B+ 

B+ 

B+ 

B+ 

B+ 

B°- 

B+ 

B+ 

B+ 

B^- 
B^- 
B^- 







P IT 

P^T] 
UJf] 

urj' 

► p+7]' 

<■ (f)n^ 
(f)rj' 



I 
I 

II 
II 
II 
II 
II 
II 
III 
III 
III 
III 
III 
IV 
IV 
V 
V 
V 
V 
V 
V 



5.7 
21 
0.75 
0.28 
0.02 
0.01 
0.46 
0.29 
6.3 
14 
6.8 
6.3 
4.5 
0.31 
0.32 
0.001 
0.040 
0.019 
0.008 
0.006 
0.001 



6.4 
23 
0.07 
0.08 
0.02 
0.001 
0.05 
0.02 
3.9 
13 
4.2 
5.5 
4.0 
0.38 
0.40 
0.0005 
0.005 
0.002 
0.0009 
0.0007 
0.0004 



7.8 
28 
1.4 
0.10 
0.06 
0.03 
0.63 
0.46 
0.89 
11 
1.0 
4.2 
3.0 
0.55 
0.57 
0.002 
0.36 
0.17 
0.073 
0.053 
0.002 



the s-quark mass decreases, as already noted in p7|j28|] . Thus, for (/i = 2.5 GeV) 
= 100 MeV, and F^^^ = 0.15, there is no problem to accommodate the CLEO data 



within the measured ±la range. As already discussed at length in refs. [p7| , p8| , these 
decay modes are dominated by the QCD penguin, and while the contributions of the 
anomaly terms are included in the rate estimates, their role numerically is subleading. 
The decay modes B^ K^rj' and B^ K^rj' show some preference for smaller values 
of ^, though this is correlated with other input parameters and at this stage one can not 
draw completely quantitative conclusions. Summarizing the B — > PP decays, we stress 
that the factorization-based estimates described here are consistent with the measured 
decay modes. All other estimated branching ratios are consistently below their present 
experimental limits. However, we do expect the modes B^ tt+tt^, B^ tt+tt'^, and 
B~^ —>■ K~^7T^ to be measured soon. 

The two observed B PV decays, B^ —>■ ujK^ and B^ —>■ uh^, = n^,K~^, show 
strong A^c-dependence as anticipated. The decay B^ un^, a class-Ill decay, has not 
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Figure 1: Branching ratio for the decays n^h^ = 7r+, K~^) as a function of for three 

different sets of form factors: BSW Model (dashed-dotted curve), Lattice-QCD/QCD-sum rules 
with central values in Table 4 (dashed curve), with the values F(^{^'^ = 0.36 and i^o^p^ ~ O-^l 
(dotted curve). The horizontal solid lines are the ilex measurements from experiment 
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Figure 2: Branching ratio for —>■ K^t]' as a function of = l/N^. The dash-dotted and 
dashed curves correspond to the choice (0) = F^'' (0) = 0.15, m,(/i = 2.5 GeV) = 100 

MeV, and Ff~"^'{0) = Fo^'^'iO) = 0.135, m,(/i = 2.5 GeV) = 122 MeV, respectively. The 
horizontal solid lines are the ±la measurements from experiment |jl|]. 
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Table 10: B —>■ PV Branching Ratios (in units of 10^^) involving b s (or \AS\ = 1) 
transitions using the BSW [Lattice QCD/QCD sum rule] form factors, with A;^ = rnl/2, p = 
0.12, 7] = 0.34, and Nc = 2,3, cx) in the factorization approach. The last column contains the 
measured branching ratio and upper limits (90% C.L.) 
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n ^ifi fn fi4l 

U.OU [U.U^I 
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ID 
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t .o ly.oj 


\ U 1 
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J\ 71 
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9 9 [9 f>l 


9 f'^ fil 
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i V 
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1 u.oo 1 
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yet been measured and the mode B~^ K^uj (a class-V decay) has a 3.9(T experimental 
significance. The branching ratios of B^ — > uoK^ and B^ —* ujtt^ are plotted as functions 
of ^ in Figs, ^and ^ respectively, showing the variations on other parameters (form factors 
and CKM matrix elements) as well. Taking the CLEO measurement B{B~^ uoK^) = 
(1.5lo;g i 0-2) X 10"^ on face value, this mode suggests that ^ < 0.1 or ^ > 0.6. The 
present CLEO upper limit B{B^ ivk^) < 2.3 x 10^^ is not yet restrictive enough. The 
branching ratio for the combined decay B^ —>■ u)h~^{h~^ = tt^, K~^) is shown in Fig ^ as a 
function of ^ for two values of the form factors Ff^^ and Ff"*'^ and two sets of values for 
the CKM parameters p and r]. The values of these form factors correspond to the BSW 
model and the upper limit in Table 4 to the Lattice-QCD/QCD sum rule case. Again, 
one sees that there is a tendency in the data to prefer larger values of the form factors. 
We note that both small values ^ — and ^ > 0.5 are compatible with data in this decay, 
with the theoretical branching ratio rising above 1 x 10~^. The value corresponding to 
the naive factorization, Ai'c = 3 (or ^ = 0.33) is definitely too low compared to the data 
on the two measured B — * PV decays. This is in line with earlier observations in the 
literature [g3,0,|l|]. 

No other B PV decays have been measured yet. However, an interesting upper bound 
B{B+ K+(t)) < 0.5 X 10"^ (at 90% C.L.) has been put by the CLEO collaboration 
0. This and the related decay B^ K^cf) are both penguin dominated and their decay 
rates are expected to be almost equal. The only worthwhile CKM-dependence is on the 
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Figure 4: Branching ratio for K^uj and as a function of ^. The legends are as follows: 

p = 0.30,7] = 0A2,Ff-'^ = 0.44 (dashed-dotted curve), p = 0.12, = 0.34, F/^^^ = 0.38 
(dashed curve). The horizontal solid lines are the ilex measurements from experiment 0. 
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Table 11: B ^ VV Branching Ratios (in units of IQ-^) using the BSW [Lattice QCD/QCD sum 
rule] form factors, with k"^ = ml/2, p = 0.12, rj = 0.34, and = 2,3, oo in the factorization 
approach. The last column contains upper limits (90% CL.) mostly from |jl]] except for the 
branching ratios for S° ^ p+p", 5° ^ p°p°, B+ p+p^, 5° ^ K*V° and B+ K*^ 
which are taken from the PDG tables 
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Wolfenstein parameter A (hence weak). However, being class-V decays, their branching 
ratios depend strongly on ^, with both having their lowest values at ^ = 0. The branching 
ratio BiB^ —>■ K^cj)) is shown as a function of ^ in Fig. |^ for A = 0.81 (dashed curve) and 
A = 0.75 (dashed-dotted curve) and the CLEO 90% C.L. upper bound is also indicated. 
This shows that values ^ > 0.4 are disfavored by the present data. In fact, taken the data 
on their face value the measured branching ratios for the decays B^ ujh^{h^ = vr"*", K"^) 
and B^ ujK^, as well as the upper bounds on the branching ratios for B^ K^cj) 
and B'^ ujtt~^ can be accommodated for a value of ^, close to ^ = 0. All other decay 
modes in Tables 9 and 10 (for the B PV case) are consistent with their respective 
upper limits. However, we do expect that the decay modes B^ — ^ P^V> ~^ P^f]' i 
B^ K*^7T^, B^ K*^7T^ and B^ p^uj should be observed in the next round of 
experiments at CLEO and at B factories. 

• There is one B VV decay mode B (pK*, for which some experimental evidence 
exists, and an averag ed branching ratio B{B (pK*) = (l.ltol ± 0.2) x 10"^ has been 
posted by the CLEO collaboration f^. The decay modes B^ —>■ (pK*^ and B^ —>■ (j)K*^ 
are dominated by penguins and are expected to be almost equal (see Table 11). They also 
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Figure 5: Branching ratio for ti^uj as a function of ^. The legends are as follows: 

p = 0.30,7] = 0.42, Fj^^^ = 0.38 (dashed-dotted curve), p = 0.12, r/ = 0.34, F^^^^ = 0.34 
(dashed curve). The horizontal solid line is the 90% CL. upper limit from experiment 0. 



belong to class- V decays, showing very strong ^-dependence (almost a factor 35!), with 
the branching ratios having their smallest values at ^ = 0. A comparison of data and 
factorization-based estimates is shown in Fig. ||. In this case, data favors 0.4 < ^ < 0.6, 
apparently different from the values of ^ suggested by the B PV decays discussed 
earlier. In fact, the branching ratios of the decays B^ — > (pK^ , B^ (j)K*^ 
and B^ (pK*^ are almost equal in the factorization approach and they all belong to 
class-V. Hence, their measurements will be rather crucial in testing this framework. 

Based on the present measurements of the B —>■ PV and B — > VV decay modes, we 
summarize that all of them belong to the class-V (and one to class-Ill) decays, for which 
the factorization-based estimates show strong sensitivity to ^. This implies that they 
are harder to predict. The classification given above, however, does not imply that the 
class-V decays are necessarily small. In fact for Nc = 2, the measured class-IV decays and 
a number of class-V B PV and B VV decays such as the ones mentioned above 
are comparable in rates (within a factor 2). For the class-V decays, the amplitudes can 
become very small in some range of ^, implying large non-perturbative renormalizations 
which are harder to quantify in this framework. Also, many class-V penguin decays may 
have significant contributions from annihilation and/or FSI, as the factorization-based 
amplitudes, depending on ^, may not dominate the decay rates. This is generally not 
foreseen for the class-I (tree-dominated) and class-IV (penguin-dominated) decays and 
most of the class-Ill decays. Hence, these decays can be predicted with greater certainty. 

Concerning comparison of our results with the earlier ones in [^,^, we note that we 
have made use of the theoretical work presented in these papers. We reproduce all the 
numerical results for the same values of the input parameters. Our decay amphtudes 
agree with the ones presented in [^, though our estimates of the matrix elements of 
pseudoscalar densities {0\u'j5u\r]^'^) and {0\d'j5d\r]^'^) differ from the ones used in pTI]. Our 
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Figure 6: Branching ratio for —>■ h^uj as a function of ^. The legends are as follows: 
p = 0.30, r/ = 0.42, Ff^^ = 0.38, Ff"*-^ = 0.44 (dashed-triple dotted curve), p = 0.30, r] = 
0.42, Ff-^^ = 0.33, Ff-^^ = 0.38 (dotted curve), p = 0.12, r/ = 0.34, Ff""'^ = 0.38, Ff"^^ = 
0.44 (dashed-dotted curve), p = 0.12, t] = 0.34, Ff^'^ = 0.33, Ff"^^ = 0.38 (dashed curve). 
The horizontal solid lines are the ilo" measurements from experiment pi]. 




Figure 7: Branching ratio for —>■ K^cj) as a function of ^. The legends are as follows: Upper 
curve: Wolfenstein parameter A = 0.81, Ff~*^ = 0.38. Lower curve: Wolfenstein parameter 
A = 0.75, Ff^^ = 0.31. The horizontal solid line is the 90% C.L. upper limit from experiment 

m. 
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Figure 8: Branching ratio for B K*(f) as a function of ^, after averaging over the and 
B^ decay rates. The legends are as follows: Upper curve: Wolfenstein parameter A = 0.81. 
Lower curve: Wolfenstein parameter A = 0.75. For the form factors, we use the BSW model. 
The horizontal solid lines represent the CLEO measurement with itlo" errors. 



expressions are given explicitly in Appendix A. The disagreement in the decay rates for 
B^ — i> p^rj and B^ p^rj' between our results and the ones given in has now been 
resolved 0. However, we do not subscribe to the notion that N(.{y + A) induced by the 
{y — A){y + A) penguin operators is different from the N^Y — A) arising from the (V^ — 
A){y — A) operators, advocated in [^] and continue to use the same N^. irrespective of the 
chiral structure of the four-quark operators. We have discussed at length the difficulties in 
predicting class- V decays some of which, in our opinion, may require annihilation and/or 
FSI effects. 

Comparison of our numerical results in the branching ratios for the B —* PV modes 
with the ones presented in requires a more detailed comment. First of all, our 
input parameters are significantly different from those of |2^. For the same values of 
input parameters, our results in charged B^ —>■ {PV)'^ decays are in reasonable accord. 
However, significant differences exist in the neutral B^ —>■ {PV)^ decay rates, which persist 
also if we adopt the input values used in ||29| . In particular, in this case we find for Nc = oo: 
B{B° p^T]) = 2.7 X 10"^ compared to 6.7 x lO'^ B{B° p^i) = 1.2 x 10"^ 
compared to 3.6 x 10^^ ||2|], B{B^ cur]) = 6.9 x 10"^ compared to 7.1 x 10~^ |29|, 
and B{B° ur]') = 1.3 x 10"^ compared to 3.6 x 10"^ [EO]. For our input values, the 



differences in branching ratios are even more drastic, as can be seen by comparing our 



results with the ones in for these decays. We have given sufficient details in our paper 
to enable a comparison of the formulae, including matrix elements of the pseudoscalar 
densities, and hence it should not be too difficult to figure out the source of the present 
discrepancy. Such details are not given in 
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• Within the present framework, we have calculated the relative importance of electroweak 

penguins in all the B PP, B — > PV and B VV decays studied in this paper. The 
decay modes where the electroweak penguins may make a significant contribution are 
shown in Table 12 where we give the ratio 

_ ^ /ii/i2)(with 07, ...,aio = 0) 

= s{B h,h,) ■ ^^^^ 

In the B — > PP case, there are five such decays whose rates show moderate dependence on 
the electroweak penguins. The decay B^ — > tt^tt^ receives significant electroweak penguin 
contribution for Nc = 3. In the class-IV B PP decays, three decays, namely -B" K'^ir^, 
B^ — > K'^T] and B^ — > K~^r] (all having branching ratios in excess of 10~^) have significant 
electroweak contributions. The presence of electroweak penguins in these decays reduces the 
decay rate by about ~ 20% to ~ 40%. 

In the B PV decays, the three class-II decays which may have significant electroweak 
penguin amplitudes are 5° p'^vr^ and B^ p'^r]^'\ Most striking among the class-lV decays 
is — > p^K^, which is completely dominated by the electroweak penguins for all values of Nc- 
This decay is estimated to have a branching ratio of 0(10"^). Measurement of this decay mode 
will enable us to determine the largest electroweak-penguin coefficient Og. In the B — > VV 
decays, the class-II decay B^ p^p^ is sensitive to the electroweak penguins. Likewise, the 
two class-IV decays, B^ p^K*^ and B^ p^K*^ are sensitive to electroweak penguins. All 
of them are expected to have branching ratios of 0(10^^) or larger, and can in principle all be 
used to determine the coefficients of the electroweak penguins. Once again, a large number of 
class-V decays show extreme sensitivity to the electroweak penguins, as can be seen in Table 
12. 



7 Stringent tests of the factorization approach and de- 
termination of form factors 

In the preceding section, we have compared available data with estimates based on the fac- 
torization approach and have already commented on the tendency of data to favor somewhat 
higher values of the form factors F^{*^, than, for example, the central values given in Table 4. 
However, as the decay rates depend on a number of parameters and the various parametric de- 
pendences are correlated, it is worthwhile, in our opinion, to measure some ratios of branching 
ratios in which many of the parameters endemic to the factorization framework cancel. In line 
with this, we propose three different types of ratios which can be helpful in a quantitative test 
of the present framework: 

• Ratios which do not depend on the effective coefficients Oj, and which will allow to deter- 
mine the form factors more precisely in the factorization framework. 

• Ratios which depend on the parameters Cj, and whose measurements will determine these 
effective coefficients. 

• Ratios whose measurements will impact on the CKM phenomenology, i.e., they will help 
determine the CKM parameters p and rj (equivalently sinct, sin/? and sin 7). 
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Table 12: Ratios of branching ratios R^^J\J defined in eq. ( ^2]) for Nc = 2, 3, oo for the form factors 
in the BSW model [Lattice-QCD/QCD sum rule method]. The horizontal lines demarcate the 
decays B ^ PP, B ^ PV and B VV. 
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7.1 Ratios of branching ratios independent of the coefficients 

We start with the ratios of branching ratios in which the effective coefficients ai, aio canceL In 
the present approach, these ratios depend on the form factors and hadronic couphng constants. 
Their measurements will allow us to discriminate among models, determine some of the hadronic 
quantities and improve the quality of theoretical predictions for a large number of other decays 
where these hadronic quantities enter. 

In what follows, we shall illustrate this by giving complete expressions for the relative decay 
widths of the decay modes in question. These expressions can be derived in a straightforward 
way from the matrix elements given in the Appendices. Then, we shall present simple formulae, 
which are approximate but instructive, and highlight the particular form factors which play 
dominant roles in these decays. Finally, we shall compare the numerical results for these ratios 
obtained from the complete expressions, which have been used in calculating the entries in 
Tables 8 - 11, and the corresponding ones obtained from the simple formulae to judge the 
quality of the approximation in each case. As practically an almost endless number of ratios 
can be formed from the seventy six branching ratios given in Tables 8-11, some thought has 
gone into selecting the eleven ratios which we discuss below. Our criterion is based on the 
theoretical simplicity and experimental feasibility of these ratios. To be specific, these ratios 
involve those decays whose branching ratios are expected to be 0(10^^) or higher, with the 
ratios of branching ratios of order one so that a reasonable experimental accuracy could be 
achieved, and whose decay widths are dominated by a single form factor. 

We start with the discussion of decay modes involving the final states tttt, pvr and pp. These 
ratios are listed below: 



Pi 



B{B^ p+p-) 



(43) 



x'f{7i,pr\Fi 



B- 



[m 



2\\2 



1 / 3x-i 



+ l)(l + a;)Mf + 



a-x^)A,A, 



where x = mp/rriB- The form factors Ai, A2 and V involve the B —>■ p transition. The function 
f{a,b) is the momentum fraction carried by the final particles, f{a,b) < 1/2. 



fia,b) 



2m| 



Since /(vr, p) ~ /(p, p) ^ 1/2— x^, and in almost all models one expects Ai ~ ^2, the expression 
given in eq. (^) gets considerably simplified. Neglecting the terms proportional to x'^ in the 
denominator, one has: 



Pi 



2\\2 



(44) 



;i + x)|Ar^(m2)p 

which is essentially determined by the ratios of the form factors Ff^'" and Ai^''. We show 
the values of the ratio Pi in Table 13 for the BSW model and the lattice-QCD/QCD sum rules 
method for both the full widths and following from the approximate relation given in eq. (0). 
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Table 13: Values of Pj's calculated with the form factors from the BSW model and the hybrid 
lattice-QCD/QCD-sum rule method. The numbers in square brackets are calculated using the 
approximate formulae derived in the text. 
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[0.47] 


0.55 [0.61] 


Pe 


1.11 


[1.21] 


1.19 [1.55] 


Pi 


1.11 


[1.21] 


1.19 [1.55] 


Ps 


1.08 


[1.14] 


0.99 [1.18] 


P9 


1.09 


[1.14] 


0.99 [1.18] 


Pio 
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There are various other relations of a similar kind. For example, neglecting the small QCD 
penguin contribution and the very small difference in phase space, we get the relations: 



B{B^ 71-71 + ) 




(45) 



As can be seen in Table 13, both eqs. (flSf) and ( [46D are excellent approximations and, for the 
two models in question, we get an almost form-factor independent prediction, namely P2 — 0.4 
and P3 ~ 0.28. It must be remarked here that one must disentangle decays from the B^ 
decays as both P2 and P3 are defined for the decays of 

In the same vein, we define the ratios P4 and P5 involving the tttt and pp modes: 

^ BiB+ ^ vr+vr") 

^- B{B+^p+p^y ^^'> 

_ B{B^ n-n+) 
= B{BO ^ p-p^) ■ 

Neglecting the QCD penguin contribution in P4 and the EW penguin in P5, which are excellent 
approximations (see Table 13), we can obtain these ratios as: 




x\l~ ml/ml) f{n,7r)\F,^-%m 



(49) 
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Neglecting higher order terms in x, we get: 



7.y \Ff-{mir 
J J {l + x)\Af-''{ml)f 



(50) 



very similar to the relation for Pi, except for the ratio of the decay constants. 
The next ratios are defined for the final states involving K*Tr and K*p. 



P7 



g(P° ^ K*+Tr-) 
~ B(PO ^ K*+p-) ' 

" B{B+ p+K*^) ' 



One can express these ratios as: 
Pe = Pr 



(51) 



(52) 



x^fiTr,K*f\Ff-^{mj,.)\ 



(l+a;)2 



(l+a;)2 



where y — mK*/'>TiB, and we have neglected the small phase space difference. Similar to the 
expression for Pi, we can derive a simple formula by dropping higher powers in x 



|Pf 



K* 



{l + x)\Ar'{ml,)\ 



(53) 



Again, neglecting the small phase space factor and the extrapolations of the form factors 
between = nip and = m^-*, the near equality Pi ~ Pe ~ P7 holds in the factorization 
assumption. 

The next ratios, called Pg and Pg, involve the final states KK* and K*K*, respectively. 
Defining 



Ps 
P9 



B{B+ K+K*^) 
B{B+ K*+K*o) ' 
-B(po ^ XOX*o) 



we now have 
PS-P9 



B{B^ ^ K*^K*^) ' 

y'\F^-^{ml.)\^\f{K,K*)/f{K*,K*)\^ 



(54) 



K7(^ + !)(!+ y)Mr? + ^^'^^ + ^ - 1(1 - 2yMrAr 

(55) 

The form factors Af, A2* , are abbreviations for Af^^* etc., and again small phase space 
differences have been neglected. Expanding in y and dropping higher order terms, we get: 



;i+z/)i^r^*«or 



(56) 
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which involves ratios of the form factors Ff^^ and Af~*^ . 

Finally, in this series we define the ratio Pio and Pu involving the states K(j) and K* 
respectively: 



Pi 



10 



- ^^^Hi^^T^. (57) 



-B(PO ^ K 

Ignoring the small phase space difference, we get 

y^\F^-^{ml)\^\f{K,ct>)/f{K%ct>)\ 



10 — -Til 



l( 3y2 



(58) 

where z = m^/mB- Expanding in y and z and again neglecting higher order terms in y and z, 
we get: 

i^lO ^ i^ll ^ 7^ M/lB^XV 2M2 • (59) 

(l+l/)l^l ("^^)P 

So, in the factorization approximation and ignoring the small extrapolation between = m\t 
and = m^, in the form factors, we have the near equality Pg ~ Pg ~ Pio — Pu. These ratios 
are all proportional to the ratios of the form factors Ff^^ and Af^^*. 

The ratios Pi, ...,Pii involve decays in which at least one of the 0~ mesons is replaced by 
the corresponding vector 1~ particle. If these particles in the decay B hih2 were heavy, 
such as D, D* , Ds, D*, one could use the large energy (1/P) expansion to derive the ratios 
Pj. We have not investigated this point and hence can not claim that these ratios are at the 
same theoretical footing as the corresponding relations involving the decays B — D[D*)7r{p), 
studied, for example, in []TB[. However, as the energy released in P ^ /11/12 decays is large. 



and no fine tuning among the various amplitudes is involved, which is the case in class-V 
decays, we think that the above relations are likely to hold. The ratios of branching ratios 
are also independent of the CKM matrix elements, therefore they constitute good test of the 
factorization hypotheses. In Table 13, we have presented the numerical values of the ratios Pj, 
i = 1, 11. This table shows that almost all the ratios are remarkably close for the two models 
used for the form factors. This, however, reflects our choice of the specific values of the form 
factors, which is influenced by the present CLEO data. In general, the ratios Pj are measures 
of the ratios of the form factors, which could vary quite significantly from model to model, and 
hence they can be used to distinguish between them. It can also be seen that in most cases, the 
simple formulae are good approximations and would enable us to draw quantitative conclusions 
about the ratios of dominant form factors in these decays. 



7.2 Determination of the effective coefficients 

In this section, we aim at measuring the effective coefficients of the factorization framework. 
To that end, we shall study some ratios of branching ratios which are largely free of hadronic 
form factors and decay constants. In general, these ratios depend on the effective coefficients Oj 
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and the CKM parameters in a rather entangled fashion. To disentangle this and gain some in- 
sight, we will have to make approximations, whose accuracy, however, we specify quantitatively 
within the present framework. 



7.2.1 Determination of the tree coefficients ai and 02 



We start with a discussion of the decays —>■ tt+tt^ and n^n^, which are on the verge 

of measurements []T|. Neglecting the electroweak contributions, which we have checked is a 
very good approximation in these decays, we can derive from eqs. ([78|) and (|80|) the following 
relation: 



^1 



TT ' TT 



TbO 



2B{B+ tt+ttO) tb^ 



ai 



ai + 02 



-2- 



ai 



ai + 02 



zi cos a cos Si + zl 



(60) 



where 



O4 + CLqRi 



ai + 02 

Here, the quantities r^o and tb+ are the lifetimes of the and B~^ mesons, which, within 
present experimental accuracy, are equal to each other. The implicit dependence on the CKM 
matrix elements in the quantity a^ + a^Ri is not very marked (see section 2). The explicit CKM 
factor is bounded from the unitarity fits in the range (at 95% C.L.): 1.4 < | V^bV^dl/I^^KTiil < 
4.6. Varying then from = 2 to Nc = 00, we get 0.08 < Zi < 0.50. This would suggest 
that one might be able to determine the quantity cos a from this ratio. However, the value of 
Zl is strongly correlated with that of the product yi = cos 5i cos a, as shown in Fig. |^ where 
the dependence of this product is shown as a function of zi, indicating the allowed range of 
Zl for assumed values of the ratio 5*1. As a result of this correlation, which is specific to the 
factorization approach, the ratio zi cos 61 cos a remains small in the entire allowed parameter 
space. The quantity zi cos^i cos a is bounded from above to lie below 0.16, which corresponds 
to using Nc = 2 and = 0.06. This is then bad news for determining the quantity 

cos a from the ratio Si but good news as far as the determination of the effective coefficients 
ai/(ai +02) from Si is concerned. Taking this as a generic case for other decays of interest, our 
best bet in the determination of the effective coefficients is to find ratios of branching ratios in 
which the quantity Zi cos 6i cos 0j (here = a, P or 7) as well as zf are both small. Within the 
factorization framework, and using the present constraints on the CKM parameters, this can be 
systematically studied. With this in mind, we shall present a number of approximate formulae 
for the ratios Si, which are expected to hold in the limit: Zi cos 6i cos (pi -C 1 and zf <C 1. 
To quantify the quality of our approximation, we shall make detailed numerical comparisons 
between the numerical results for Si, obtained with the complete expressions for the respective 
decay widths, and the ones following from our approximate formulae. 

There are some ratios of branching ratios in which, within our theoretical framework, the 
factors Zi cos Si cos (pi are large, or else the CKM dependence of the ratios factorizes in a simple 
way. We shall use these ratios to determine the CKM parameters in non-leptonic two-body 
decays B hih2- This kind of analysis has already been suggested in the literature ]45|J27| , |5D[] . 
We add a number of interesting decay modes to the cases already studied in the literature and 
make quantitative predictions for them in the present model. 
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Figure 9: yi — cos^icoso; as a function of Zi in the factorization approach. The dotted, 

dashed-dotted and dashed curves correspond to iV^ = oo and iKib/V^tl = 0.11, = 3 and 
iKb/Kbl = 0.08 and Nc = 2 and iKfc/Kbl = 0.06, yielding in the BSW model the values 
5*1 = 2.07, Si = 0.94 and Si = 0.59, respectively. The two vertical lines indicate the bounds on 
zi from our model and the CKM unitarity fits 0.08 < zi < 0.50. 



Returning to the determination of the coefficients a^, we note that a ratio similar to Si can 
be defined with the pp final states: 



^2 = 



P P~ 



Tbo 



2B{B+ ^ p+pO) TB+ 



ttl 



where 



Z2 



Oi + 02 



Oi + 02 



-Z2 COS a COS ^2 + z. 



(61) 



04 



Ol + 02 



This is not expected to exceed its maximum value 2;^"^ = 0.26, the central value being around 
Z2 — 0.08. Hence, one could use an approximate formulae for and 5*2 by keeping the dominant 
term arising from the tree contributions (setting tbo — Tb+)'- 



^1 = 
^2 = 



2B{B+ 
B(B^ 



> tt+ttO) 

pV) 



Oi + O2 
Oi 



(62) 



(63) 



2B(B+^p+pO) V01 + 02. 
Likewise, neglecting the penguin contributions, which give only several percent uncertainties, 
the value 02/01 can also be measured from the following ratios. 



^3 

^4 = 



^ 2B{B+ p+7r°) 
~ B(S0^p+7r-) 

_ 2B{B+ ^ TT+p") 
B{B^ TT+p-) 



1 + 



102 y 

X aiJ 



~ 1 + a; 



02 

Oi 



(64) 
(65) 



where the quantity x = (f pFf (f^Ag'''') can be measured by measuring the ratio P3. 
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Table 14: The ratios Si calculated using the indicated values of N^. and different values of p 
and 77. The values are calculated using the approximate formula (Approx.) derived in the text 
also. 









Nc = 


2 


Nc = 3 


Nc = oo 






0.06 


0.08 


0.11 


0.06 


0.08 


0.11 


0.06 


0.08 


0.11 


Si 


Exact 


0.59 


0.66 


0.68 


0.83 


0.94 


0.95 


1.81 


2.03 


2.07 




Approx. 


0.64 


0.64 


0.64 


0.91 


0.91 


0.91 


1.97 


1.97 


1.97 


S2 


Exact 


0.60 


0.64 


0.65 


0.85 


0.90 


0.91 


1.84 


1.95 


1.98 




Approx. 


0.64 


0.64 


0.64 


0.91 


0.91 


0.91 


1.97 


1.97 


1.97 


^3 


Exact 


1.33 


1.32 


1.32 


1.09 


1.09 


1.09 


0.74 


0.75 


0.75 




Approx. 


1.29 


1.29 


1.29 


1.06 


1.06 


1.06 


0.71 


0.71 


0.71 


5*4 


Exact 


2.41 


2.20 


2.13 


1.40 


1.23 


1.17 


0.32 


0.23 


0.19 




Approx. 


2.13 


2.13 


2.13 


1.20 


1.20 


1.20 


0.22 


0.22 


0.22 


^5 


Exact 


0.55 


0.97 


1.96 


0.37 


0.66 


1.33 


0.16 


0.28 


0.56 




Approx. 


0.55 


0.95 


1.89 


0.38 


0.66 


1.31 


0.17 


0.29 


0.58 


Sq 


Exact 


3.07 


5.46 


11.01 


1.95 


3.47 


7.00 


0.75 


1.34 


2.71 




Approx. 


3.10 


5.30 


10.56 


2.06 


3.53 


7.03 


0.86 


1.46 


2.92 


S7 


Exact 


1.97 


3.73 


7.62 


1.77 


3.35 


6.84 


1.49 


2.82 


5.76 




Approx. 


1.99 


3.40 


6.78 


1.87 


3.19 


6.36 


1.68 


2.88 


5.74 


Ss 


Exact 


1.84 


3.47 


7.10 


1.65 


3.12 


6.37 


1.39 


2.62 


5.36 




Approx. 


1.99 


3.40 


6.78 


1.87 


3.19 


6.36 


1.68 


2.88 


5.74 


Sg 


Exact 


0.32 


0.65 


1.33 


0.31 


0.62 


1.27 


0.28 


0.57 


1.17 




Approx. 


0.36 


0.61 


1.22 


0.35 


0.59 


1.18 


0.33 


0.57 


1.13 


5*10 


Exact 


0.22 


0.18 


0.13 


0.15 


0.14 


0.13 


0.09 


0.12 


0.18 




Approx. 


0.14 


0.14 


0.14 


0.13 


0.13 


0.13 


0.11 


0.11 


0.11 


5*11 


Exact 


0.37 


0.17 


0.06 


0.28 


0.15 


0.07 


0.20 


0.16 


0.12 




Approx. 


0.26 


0.15 


0.07 


0.25 


0.14 


0.07 


0.23 


0.13 


0.07 



7.2.2 Determining the penguin coefficients 

Concerning the coefficients 03, Oq, we recaU that the dominant QCD penguin amplitudes are 
proportional to 04 and Og. The others (03 and a^) enter either as small corrections in class-IV 
decays, or else enter in class-V decays, which in most cases are highly unstable due to large 
cancellations in the respective amplitudes, hence rendering this exercise not very trustworthy 
for determining the smaller coefficients. In view of this we concentrate on relations involving 
the QCD-penguin coefficients and a^. For this purpose, quite a few class-lV decays listed in 
Tables 8-11 suggest themselves. Here, we take the ratios between some of the representative 
decays from this class and from class-I or class-Ill decays. These ratios and their approximate 
dependence on the coefficients of interest are as follows: 



S,= 



_ 2B(B+ tt+tt") 




ai + a2 



CI4 + CLqR^ 



(66) 
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Here, the quantity R; 




(67) 



(68) 



(69) 



(70) 



As is obvious from 



the formulae given above, the determination of the effective coefficients through these ratios is 
correlated with the values of the CKM factors, which in all cases are given essentially by the 
ratio iKb/V^sl — iKb/Kbl - 0.08 ± 0.02. We expect that the CKM matrix element |Kb/Vd.| 
will be very precisely measured in forthcoming experiments. Hence, a better use of these ratios 
is to determine the effective coefficients. To give a quantitative content to the approximations 
made in reaching the simple expressions for Si, i = 1, 9, we display in Table f4 the numerical 
values of these ratios, together with the ratios 5*10 and 5*11 discussed below, as a function of 
l^b/Kbl) taking a rather generous error on this quantity in the range 0.06 < iKib/Kbl < 0.11, 
for three values of A^c- The rows labeled as "Exact" are the results obtained by using the 
complete amplitudes and the rows labeled as "Approx." are the results following from the 
simple relations given above for these ratios. As one can see, these formulae are quite accurate 
over a large parameter space, with the deviations mostly remaining well within 10%. One can 
also check that the ratios - Sg for the complete result scale almost quadratically with Kb/ Kb, 
as follows from the simple formulae, which shows that the CKM dependence displayed in the 
approximate formulae is actually quite accurate. 

Concerning the measurements of the electroweak coefficients, a7,...,aio, we recall that the 
dominant contribution of the electroweak penguin amplitudes is proportional to ag. The rest of 
the electroweak coefficients are either small or they enter in combinations which render them 
very sensitive to the variation in N^. It is instructive to consult Table 12, where the decays 
in which electroweak penguins may make a significant contribution to the branching ratios are 
listed. In line with our argument, we will concentrate only on class-IV penguin-decays, and 
pick up the decay mode — > p'^K'^ as an illustrative example. To that end, we define the 
following two ratios involving a class-I and a class-IV processes, dominated by the tree and 
QCD-penguins, respectively.: 



'10 



B{B+ - 
2B{B^ 



> p^K^) 



B{B^ p-n+) 



fpFl 



2 






a4 



^^^(m2) 



'''{ml) 



Og 

ai 



(71) 



We show the numerical values of these ratios in Table 14 for the three indicated values of 
the ratio iKtt/V^bl, both for the exact and approximate cases. The approximate relations are 
reliable over most part of the parameter space. Other similar ratios can be written down in 
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Figure 10: y\2 = cos5i2Cos7 as a function of 2:12 in the factorization approach. The dotted, 
dashed-dotted and dashed curves correspond to Ai'c = 00 and |V^tif,/Kfe| = 0.11, A^^^ = 3 and 
\Vub/Vcb\ = 0.08, and Nc = 2 and \Vub/Vcb\ = 0.06, yielding in the BSW model the values 
5*12 = 0.46, 5*12 = 0.91 and 5*12 = 1.12, respectively. The two vertical lines indicate the bounds 
on 2:12 from our model and the CKM factors discussed in the text, yielding 0.15 < 2:12 < 0.29. 



a straightforward way. Measurements of the ratios Si - Sn will overconstrain the coefficients 
04, ag and og, testing both the factorization hypothesis and determining these crucial penguin 
coefficients. Note that Sio depends only slightly on the CKM factors, compared to the others 
discussed above, and Si to S'4 do not depend on |Ki;,/Kb| when we use the approximations in 

eqs. m-m- 



7.3 Potential impact of B ^ hih2 decays on CKM phenomenology 

(i) B nK channels: 

In this subsection, we consider the ratios of branching ratios which can be gainfully used 
to get information on the CKM parameters. The most celebrated one in this class is the 
ratio discussed by Fleischer and Mannel recently [^, involving the decays —>■ K~^7t~ and 
B~^ — > K^TT^. Ignoring the electroweak penguin contribution, which is estimated to be small 
in our model, one can write this ratio as: 



with 



S 



B{B^ 



12 



2^12 



TT 



\T\ _ 



1 — 2^:12 cos 5i2 COS 7 + z 



2 

125 



(72) 



VtbVtl 



ai 



The branching ratios involved in S12 have been measured by the CLEO collaboration and 
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their values can be seen in Table 8. The ratio S12 itself has the following value: 



S 



12 



0.65 ±0.39 



(73) 



For the central values of the CKM parameter (p = 0.12, 77 = 0.34), the value of S12 is found 
to be 0.80 < 5*12 < 1.0 varying Nc and using the two form factor models displayed in Table 8. 
However, varying the CKM parameters in their presently allowed range, we find 0.46 < 5*12 < 
1.12, where the lower and upper values correspond to |\4f,/K;fe| = 0.11 and |\4fe/V^f,| = 0.06, 
respectively. The ratio 5*12 is, formally speaking, very similar to the one defined for the ratio 
5*1. However, the difference between 5*1 and 5*12 is that the product 2:12 cos5i2 cos 7, as opposed 
to the corresponding quantity Zi cos^i cos a in 5*1, is not small in the allowed region of Z12. The 
range 0.15 < Zu < 0.29 is estimated in the factorization approach varying the CKM matrix 
element ratio in the range 0.013 < IK^V^JJ/I^tb^tsl < 0-023 and Nc. This is shown in Fig. 0. 
Hence, the ratio 5*12 and its kinds, discussed below, do provide, in principle, a constraint on 
cos 7. This figure also shows that the ratio S12 is in quite good agreement with the measured 
ratio by CLEO. 



In the context of the factorization models, the CLEO data was analyzed in [g^ and it was 
shown that theoretical estimates in this framework are in agreement with data. The ratio S12 
(called Ri in [^) provides a constraint on the CKM parameter p (equivalently cos 7). Taking 
data at ±lcr value, the CLEO data disfavored the negative-p region. The allowed values of this 
parameter resulting from the measurement of 5*12 were found to be in comfortable agreement 
with the ones allowed by the CKM unitarity fits. In addition, the dependence of 5*12 on the 
CKM parameter rj was found to be weak. This overlap in the value of p following from the 
analysis of the ratio Su in the factorization approach and from the CKM unitarity fits has 
also been confirmed recently in We show here the ratio S12 plotted as a function of 

cos 7 for A^c = 2,3 and 00 and fixed value of the ratio \Vub/Vcb\ = 0.08 in Fig. |ll]. The form 



factor dependence of this ratio is rather weak (as can be seen in Table 8) and for the sake of 
definiteness we display the result for the BSW form factors. It is seen that for all values of A^^^, 
the CLEO data provides a constraint on cos 7, which is compatible with the one allowed by the 



CKM fits, yielding 32° < 7 < 122° [0. This is in line with what has already been reported in 
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The ratio S12 given in eq. ( [7^ ) is a generic example of the kind of relations that one can get 
from the ratios of branching ratios in which the quantity Zi cos 6i cos 7 is not small. We have 
argued, in line with |^, that the factorization model gives an adequate account of S12. We 
discuss below some related ratios, which, once measured, could be used to determine cos 7 as 
well as further test the consistency of the factorization approach. 

(ii) Ratios for B nK* modes: 

One can define analogous to eq. (|72|), the ratio S'13, involving the decays 5° n^K*^ and 



TT 



+ . 



'S'13 = 



i3(50 



TT 



B{B- 



with 



TT" 



^13 



-K 



|T| 



1 — 2^13 cos (5i3 COS 7 + 2:^3, 



(74) 



ub 



V* 

' us 



VtbVt*s 

Using 0.013 < \VubV*^\/\VtbV*^\ < 0.023, and from A^", = 2 to A^'^ = 00, we get 0.30 < z^^ < 0.60, 
indicated in Fig.lKl The ratio S'13 is plotted in Fig. O as a function of cos 7 for three different 



ai 
2 to Nr 
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Figure 11: 5'i2 as a function of cos 7 in the factorization approach. The dotted, dashed-dotted 

and dashed curves correspond to = 00, = 3 and = 2, respectively. The horizontal lines 
are the CLEO (ilo") measurements of 6*12. The two vertical lines correspond to 32° < 7 < 122°. 
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Figure 12: 1/13 = cos ^13 cos 7 as a function of 2:13 in the factorization approach. The dotted, 
dashed-dotted and dashed curves correspond to Nc — 00 and |T4fe/K;6| — 0.11, Nc — 3 and 
iKfc/V^cbl = 0.08, and Nc = 2 and iKb/K&l = 0.06, yielding in the BSW model the values 
5*13 = 0.49, 5*13 = 0.95 and 5*13 = 1.37, respectively. The two vertical lines indicate the bounds 
on Zi3 from our model and the CKM factors 0.30 < 2:13 < 0.60. 



values of Nc and iKife/Kfel- When measured, this ratio will provide a constraint on the phase 
cos 7. Varying the CKM parameters and Nc in the indicated range, we find the ratio Sis to lie 
in the range 0.49 < Sis < 1-37. The upper bound is larger than the one for S'12 given earher, 
reflecting that the QCD-penguin contributions in the two ratios are similar but not identical, 
(iii) Ratios for B pK modes 
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Fi gure 13: Si^ — as a function of cos 7. The dotted, dashed-dotted and dashed hnes 
correspond to results with = 00, Nc = 3 and = 2, respectively. The two vertical hnes 
correspond to 32° < 7 < 122°. 



S 



14 



^ p-K+) 



1 — 2^14 cos 5i4 COS ■y + z 



2 

145 



(75) 



with 



\T\ 



ai 



VtbVtl a4 + 06^4 

The central value of the quantity ^14 IS Zi4 ~ 5.07. However, being very large, the ratio 5*14 
implies that the branching ratio in the denominator is appreciably smaller and perhaps lot 
more difficult to measure. In view of this, we are less sure of its utility of the ratio Sm in the 
foreseeable future. 

(iv) Ratios for B pK* modes 

Finally, we note that the ratio S'15 defined below provides, within our model, a very similar 
constraint on cos 7 as the one following from the ratio 5*13: 



5'l5 = 



1 — 2Zi5 COS 5i5 COS 7 + 2 



(76) 



13{B+ p+K*^) 

where 2:15 = ^13 and ^15 = 613. This will be a further test of the factorization Ansatz. 

Finally, in conclusion of this section, we mention that a method of measuring the CKM 



matrix element ratio 
ref. [| 
effects- 
below: 



\Vtd/Vt 



ts\ 



using exclusive non-leptonic B decays has been proposed in 
. Some of these ratios have modest theoretical uncertainties due to SU(3) breaking 
These relations hold in the factorization framework as well, and we list a few of them 



B{B+ -> K+K^) _ B{B+ K+K 



B{B^ 



K*+K^ 



'K^) B{B+ - 
B{B+ K* 



B{B+ p+K^) 



(77) 



B{B+ ^ p+ K 



td 



ts 



46 



8 Summary and Conclusions 



We have presented estimates of the decay rates in two-body non-leptonic decays B — > hih2 
involving pseudoscalar and vector hght hadrons in which QCD and electroweak penguins play 
a significant role. This work partly overlaps with studies done earlier along these lines on 
branching ratios, in particular in [p^^ -|3T| . We make use of the theoretical framework detailed in 



[27, 28 1 but we think that this is the most comprehensive study of its kind in the factorization 
framework. 

Using the sensitivity on Nc as a criterion of theoretical stability, we have classified all the 
decays B hih2 into five different classes involving penguin and tree amplitudes. This extends 
the classification of tree amplitudes en vogue in the literature [§|,|18|. We hope that the detailed 
anatomy of the decays B hih2 presented here, in particular concerning the QCD and elec- 
troweak penguins, will serve to have a more critical view of what can be reasonably calculated in 
the factorization framework and what involves a good deal of theoretical fine tuning. Following 
the classification discussed here, we think that class-I and class-IV decays, and probably also 
class-Ill decays, can be calculated with a reasonable theoretical uncertainty, typically a factor 
2. However, most class-II and class- V decays deserve a careful theoretical reappraisal to es- 
tablish the extent of non-factorizing contributions. In particular, we have outlined the pattern 
of power suppression in annihilation contributions to two-body non-leptonic B decays. Being 
suppressed by m^/m^, the annihilation contributions are small in the decays B PP but 
since this suppression is only m\/m% 'm B ^ PV and B —>■ VV decays, in specific cases this 
can be easily overcome by the favorable effective coefficients. Hence, annihilation contributions 
can be significant in some B —>■ hih2 decays involving vector mesons. 

Our results can be summarized as follows. 

• The recently measured decay modes B^ —>■ K^n^,B^ — > K~^r]',B^ — > K^r]',B^ — >■ 
7r~^K^, and B~^ —>■ ujK^ can be explained in the factorization framework. The first four 
of these belong to the QCD-penguin dominated class-IV decays, which we argue can 
be reliably calculated. The last belongs to the iVc-unstable class-V decays, which may 
receive significant FSI and/or annihilation contributions. Taken the present theory and 
data on face value, all measured decay modes are consistently accommodated, with some 
preference for ^ = l/N^. < 0.2. Data on the combined decay modes B (pK* prefers 
somewhat higher value for ^. However, we caution against drawing too quantitative 
conclusions at this stage. 

• A number of decays is tantalizingly close to the present experimental upper limits. We 
think that with O(IO^) B/B hadrons, available in the next three to five years, a good 
fraction of the seventy six decay modes worked out here will be measured providing a 
detailed test of the factorization approach. 

• To further quantify these tests, we have put forward numerous proposals which involve 
measurements of the ratios of branching ratios. Carefully selecting the decay modes, one 
could determine the effective coefficients 01,02,04,00 and 09 from data on B ^ hih2 
decays in the future. A consistent determination of these coefficients will greatly help 
in developing a completely quantitative theory of non-leptonic B decays. Leaving out 
02 from this list, which depends significantly on Nc, we do not expect that the rest 
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will be greatly modified by non-perturbative effects. It will be difficult to quantitatively 
determine the smaller penguin coefficients not listed explicitly. 



We have proposed a number of ratios involving the decays B hih2, relating the final 
states in which a pseudoscalar meson is replaced by a vector meson. They will help in 
determining the form factors for the various decays considered here. While these relations 
are derived in the factorization approach, perhaps their validity is more general. 

The current and impending interest in two-body non-leptonic decays for the CKM phe- 
nomenology is illustrated, arguing that they provide potentially non-trivial constraints on 
the CKM parameters. While ultimately not competitive to more precise determinations 
of the CKM parameters from B^-B^ mixings and radiative and semileptonic B decays, 
they are of current phenomenological interest as the constraints following from them are 
already complementary to the ones from the CKM unitarity fits. 

Finally, within the factorization framework which gives an adequate account of the present 
data on decay rates, it will be instructive to study direct and indirect CP violation in 
all two-body non-leptonic B decays discussed here. We hope to return to this in a 
forthcoming publication . 



Acknowledgments 

We thank Christoph Greub and Jim Smith for helpful discussions on various aspects of theory 
and data discussed in this paper. We also thank Jim Smith for critically reading the manuscript. 



A Matrix elements for B decays to two pseudo-scalar 
mesons 

(1) b ^ d processes: 



X {VubV*aai-VtbV;a[a4 + aio + {aQ + as)Ri]}, (78) 



with Ri ^"^ 



{mi,-mu){mu+md) ' 



M(5°^7rV) = ^^/.Fo^--(m^)(m|-m^){K.K>2 

13 3 11 
+VtbVtd[a4 - -aw + -aj - -ag + (ag - -as)R2\ | , (79) 

with i?2 = 7 ^r^^ — 7. After squaring of the matrix element, the decay rate should be 

divided by 2, for the symmetric factor of identical particles in the final states. 

M{B-^n-n') = -^^^F,^-{ml){ml - ml) 

X ^VubV*a{ai + a2) - VtbVta X ^[ag + aio - aj + a8R2]Y (80) 
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M{B- 



TT 1] 



.Gf 



— i- 



1 1/ X / 

-aio + 2a3 - 2a5 + -(ag - ay) + (ag 



+ (a3 - as + ag - a^)-^ + [a^-a^- -(ag - ay 

/„(') 



1 



fu 

/„(') 



(') 

3 



1 - 



fu 
/^(') 



with 



„(') 



(mt-md)ms 



The definitions of the decay constants involving rj and t]' are as follows: 



(0|M7^75M|r7W(p)) = zKop, , (0|s7^75s|r7W(p)) = .fs^^ ^ {0\ci,l5c\v^'Kp)) = ^f'ioV, ■ (82) 



The quantities /^(,) and /^(,) in the two-angle mixing formalism are: 

In' 



sin ^8 + COS 6*0, 



(83) 



(84) 



We shall also need the matrix elements of the pseudoscalar densities for which we use the 
following equations: 

(0|s75s|r7) ' 
(0|s75s|r/') f^, ' 



(85) 



which differ from the corresponding equations in which have been sometimes used in the 
literature. In the approximation of setting /§ = /o, and 6*8 = 6q, the relations given above, 
however, agree with the results derived in |7^. The results for the densities (0|s75s|?7') and 



(0|s75s|?7) have been derived in which we use here: 

(0|s75s|?7') 
(0|s75s|r7) 



—I 



—I 



2ms 
2m„ 



(86) 



We point out that the anomaly contributions have been taken into account in deriving these 
expressions. They are numerically important. The relevant form factors for the B ^ rj' and 
B ^ 1] transitions are: 



pB-r,' ™ /sinA , cos^\ 



' cos sin 6*0 ' 



(87) 
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The mixing angles that we have used in the numerical calculations are 6g = —22.2°, 6q 



[71 



A^'^) = -^^^Fr''\rnl)irnl-rnt,,){V^,V:,a, 



1 / 1 NO 3, 



/«{') 



2 



j„(/) 



a4 + 2a3 - 2a5 + ^(ag - 07 - aio) + (ag - ^as)^'^ | 1 



/'^(') / 1 \ /*(/) 

+ (03 - as + ag - 07)-^ + (as - ^5 - 77(^9 - a?) I 7^ 



M(S° ^ 7]r]') 



04 + 2a3 - 2a5 + ^(og - 07 - aio) + (oe - ^a8)-R3 f 1 



+ (03 - 05 + ag - 07)-^ + (^fls - as - -(ag - a7)j -j^ | 



a4 + 2a3 - 2a5 + ^(ag - a7 - aio) + (^6 - ijO,fi)B!^ ^1 
+ (a3 - as + ag - ^r)-^ + (^«3 - 05 - - aj)^ | . 



a4 + 2a3 - 2a5 + i(ag - a7 - aio) + (ae - ]jCi^)B!^ il 



+ (% - ^5 + ag - a?)^ + Us - as - -(ag - a7) 



Jri' 



The matrix elements for M{B^ rji]) are the same with the above equation with rj' 
(2) b ^ s processes: 



M(5° 



K- 



X {KfeK^ai - V't6V^*[a4 + aio + (ae + a8)-R4]} 
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with i?4 



2ml 



(m6-m„)(m„+ms) ' 



1 1 AD ' 

04 — -aio + (ae — -^O'SjRb 



(92) 



with it^s 



2m2 _ 



(m6-md)(md+ms) ' 



-ytbVZ[ai + aio + (ae + a8)i?4]} (93) 



with i?^'^ = 



2m,' 



(') 



(') , 



-Vt6V;*[a4 + aio + (^6 + 

\/2 2 



ts 



203 - 205 + ^(ag - 07) - (ae - ^a8)-R6 



+ (as - 05 + ag - 07)7^ + 

/r,(') 

^03 - 05 + a4 + ^(07 - og - oio) + {ae - ^as)-^! 



(mb—ms){ms+ms) ' 



(94) 



M(5° ^ i?%W) 



■ Gp B^rf'K 2 \ / 2 2 \T/T/* 

i'n^K){'mB-m^('))VtbVt, 



1 1 

04 — -Oio + [0,6 — -^dsjR^ 



2a3 - 2a5 + -(ag - 07) - (oe - ■^o,8)R6 



(') 



+ (03 - 05 + ag - 07)-:^ + 

/„(') 



(0,3 - a5 + a4 + ^(ay - ag - aio) + («6 - ^«8)-R6 



fu 



(95) 
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(3) Pure penguin processes: 

M{B- ^ n-K') = i^fKFi-''{ml){ml - ml)Vt,VZ {04 - \a,o + {a, - \as)R,] . 



(96) 



with Rf 



2-m? „ 



(m6-ms)(md+m3) ' 



M(S" ^ X"i?") = i^/i^Fo^-'^(m^o)(m^ - m],,)Vt,V:^ |a4 - ^Oio + (ag - |a8)i?r| • (98) 

B Matrix elements for B decays to a vector and a pseudo- 
scalar meson 

(1) 6 — > d processes: 

M(SO ^ p-7r+) = V2GFf,Ff^-{ml)m,{e ■ p^) {KbKTd^i - ^tb^tdK + aio]} • (99) 



1 



1 



M(S0^p+7r-) = v^G'^/,<-^''(m^)m,(e-p,){K6K>i 

-^t6^td[«4 + Clio + (oe + a8)(3i]} , 



with Qi = , , 



M(5° ^ 7rV°' 



+ 



fpFf-^ml) |v;;.K>2 + T46T4*J«4 - ^aio - ^(ar + ag)]} 



04 - -aio + (^6 - ■7:a8)Q2 + ^(a? - og) 



2m^ 



with (m6+md)(md+md) 



(100) 



(101) 



M{B- ^ 7r-p°) = Gprnpie ■ p,) (UA^^'{ml) {KbOi - V^bV^t^K + aio + {a, + a8)gi]} 

+ f,Ff--^ml) {v;6K>2 - ^t6KJ-a4 + ^aio + ^(a^ + ag)]}) . (102) 



M{B- ^ p-7r°) = Gprnpie ■ p^) (uA^^%ml) \ ^6^2 - Vt^V, 



'td 



-04 + -oio 



1 3 
+ (-«6 + -afi)Q2 + -(ag - a?) 

+ UFf-^'iml) {VubV:,a, - VM[a, + aio]}) ■ 
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(103) 



G 

-Vtiyu 



11 3 ■ 

-a4 + -aio + (-as - a6)(52 + -^^iag - a?) 



a4 + 2(03 + as) + - (ay + ag - aio) 



(104) 



M(S- ^ TT-cu) = Gprn^^e ■ p^) [UA^^'^^ml) {K6V;>i - VftV^^K + aio + (^6 + a^)Qi]} 



td 



a4 + 2(a3 + as) + - (ay + ag - aio) 



(105) 



M{B- ^ p-ryW) = V2GFmp{e ■ p,,)) (jX^"^'" H) {KbKVi - ■^"t6l4d[«4 + aio]} 

/«(') 



(') 



I 



1 1 / /"{/) 

a4 + 2a3 - 2a5 + -(ag - ay - aio) + (oe - 2^8)Q? I 1 - jj- 



+(a3-a5 + ag-ay)^+ (a3-a5-;^(ag-ay)) ^ I J, (106) 
/r,(') \ z / /^(,) J J y 



where Q 



M{B'^p%^'^) = G^m,(e-p,(,))(/,Ff^^*'(mJ){K,i;>2 



1 3 , 
-a4 + -aio + 2('^9 + ay) 



- /;(oAr'(mJ')) ^ ^«''K.>2 + v;6K>2^ 



1 



a4 + 2a3 - 2a5 + -(ag - ay - aio) + («6 - ■;^«8)Q, 



2 

+(a3 - as + ag - ay)-^ + (a3 - as - -(ag - ay) ) 



(') 
3 



fu 

(107) 
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04 + 2(03 + 05) + -(07 + Og - aio) 



I /,,(') 



td 



a4 + 2a3 - 2a5 + ^(ag - ay - Oio) + (ae - ]jas)Qf ( 1 - 



fu 



+ (03 - 05 + 09 - a7)4r" + («3 - 05 - ^(Og - 07)") ^ I I . (108) 
■/■nC) V ^ / /^(,)J J J 

{2) b —>■ s processes: 

M(S° ^ K*-7r+) = V2GFfK*Ff-^-(ml.)mK*(e ■ p^) {K6K>i - ■«^t6V,:[a4 + aio]} . (109) 



V2GFfKA^^'{mj,)m,{e ■ pk) {K6K>i 
-"Vt6V^t5[«4 + oio + (oe + as)Qi]} , 



(110) 



with (54 



-2mi 



VubV*,a2 - VtbV*-{ag - ay) 



^4 — -CtlO 



(111) 



M(B°^i?V°) = G^m,(6-p;,)|/K^^^''(m2,o)Wt: 



1 1 
04 - -oio + (oe - 2"8)(35 



VubV*sa2 - VtbVtl X -(ag + 07) 



(112) 



with Qs = 



(mb+md)(md+ms) ' 



M{B- ii-^-Tr") = GprnK-ie-p^) 



UA^-^\ml) [v^bV:,a2 - V^,V;, x ^(ag - ay)] 



+ fK*F,^-*-{ml.) {K6K>i - 146l4:(a4 + aio)}] ■ 



(113) 



M{B- K-p^'' 



= Gpnipie ■ Pk) [JkAq ^^(m^) {KtK^i 

-VtbV*,[a4^ + aio + (fle + a8)<34]} 
+ f,Ff-''{ml) [VubV:,a2 - VtbY^ x ^(ag + a,) 



(114) 
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a4 



^aio + (ae - ^as)Q5 



2(03 + 05) + -(09 + 07) 



(115) 



M{B- K-(jj) 



Gprn^ie ■ Pk) /i^A^"*'^(m^) {KfeK*s«i 

-^tfe^tt(«4 + aio + (cie + a8)Q4)} 
+ UFf^^'iml) {k6C«2 - 1^*6 V,: (2(03 + as) + ^(ag + ay)) 



(116) 



M(B" 



X*-7yW) = V2GFmK'{e-pB)(^fK'Ff^'''\ml){V^bV:,ai-VtbV:,{a4 + a,o)} 



with Q^'^ = 



1 1 

2(03 - as) + -(ag - a?) - (ae - -a8)(36 



(') 



+ (aa - as + ag - ay)-^ + 



(as - as - ^(ag - aj) + a4 - ^aio + (ae - ^a8)Qe^ 



fu 



2m" 



(') 



(mi,+ms)(ms+m3) ' 



(117) 



M(S° ^ ^*^W) = v^G^m^. (e ■ pb) i -fK-^F^^'^'" {ml)Vt,V, 



ts 



a4 — 2*^10 



2(a3 - as) + ^(ag - ay) - (ae - ^a8)(5j^+ 

/«(') 



+ (ag - as + ag - ay)^ + 



^a3 - as - ^(ag - ay) + a4 - ^aio + (ae - ^a8)Q6^ 



;/)\ f^i') 



fu 

/„(') 



(3) Pure penguin processes: 

M{B- ^ 7r-i?*°) = -V^G'f/K*Ff^"(m2,,)mK*(e -p^VibK^ 



a4 — -aio 



(118) 



(119) 
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M{B- ^ p-K^) = -V2GFfKA^-f{mlo)m,{e-pK)VtbV* 



■ (120) 



td 



04 — -^aiQ 



(121) 



M{B- K*-K^) = 



X VtbV^ 



td 



04 - -aw + (oe - 2'^8)Q7 



with Q 



-2m- 



{rni,+ms){md+ms) ' 



(122) 



M(S- ^ 7r~(/)) = -V2M{B^ 



(123) 
(124) 



M(S° ^ 77W0) = -V2GFUF,^-'''"\ml)m^{e ■ p^^)Vt,V;, jag + 05 - ^(ay + 09)}. (125) 



M(s- ^ /s:-(/)) = 

M{B°^K^(f)) = 



X 



jaa + a4 + 05 - -(07 + 09 + oio)| . 



(126) 



C Matrix elements for B decays to two vector mesons 

(1) 6 — > d processes: 

M{B^^p-p+) = -i9Zf^m,{{e+-e^){mB + m,)Af-'{ml) 



-(e+ •pB)(e_ -pb)- 



{vriB + mp) 

X {KbK>l - ^i6^w[«4 + Oio]} • 



'^^lxva(3^—^-\-PBP+' 



(niB + 



(127) 
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(ei ■pB)[e2-PB)-f \ ^ - iei^^ai3ei€2PBP2 



X 



f 13 3 1 

|K6K.>2 + V^foV^^K - 2«io - - 2«9]| • (128) 



M{B- p-p^) 



ml) 



-(eo ■PB){e- -Pb) 



X 



{niB + mp) 
3, 



P' u cx p ^ P^ 

l^pua/3^-^oPBPo 



{niB + nip) 



3 ^ 

yubV*d{ai + 02) - VtbV^*a X -[ay + ag + aio]| . 



(129) 



M(S°^u;u;) = -i^/a,m^{(ei-e2)(mB + m^)Af^"(m^) 



-(ei • PB){e2 ■ Pb)-, : r - i^u.^ape'le^p^p'^ 



X 



[mB + m^) [mB + m^) J 

|K6V;>2 - VtbV;a[a4 + 2(03 + as) + ^(ar + ag - aio)]| . (130) 



MiB'^p'u) = -i^f^mp{{eo-e^){mB + m^)A^-^{ml) 

2Af^-(m2) ^ ^2V^^^"(m2^ 

13 1 

K6V;>2 - VtbV;a[-a4 + -aio + 2(^7 + a9)]| 

i^^Um^ {(eo • eu;)(mB + mp)Af'^^(ml) 
-(eo • PB){e^ ■ Pb) ,^ , , - ^e^.a/3e^ie^p>^ ^ 



X 



+ i- 



X 



{niB + mp) (m^ + mp) J 

|K6Kr(i«2 - VtfeV^^K + 2(a3 + as) + ^[ay + ag - aio]| . (131) 



G 

M{B-^p-u) = -i^/pmp{(eo-e_)(mB + m^)Af^^(m2) 

, w . 2>lf-(m^) . , _ , 2y^-(m^) 
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X {K(,K>i - ytbV:a[aA + aio]} 

- i-ffujm^{{(^o-(^-){rnB + mp)Af^^{ml) 



X 



(niB + nip) 



i^HuaP^O^-PBP- 



{niB + nip) 



04 + 2(03 + 05) + -(07 + ag - Oio) 



(2) 6 — > s processes: 



-(e+ •pB)(e_ -ps) 



2A 



m 



K* 



{iriB + nip) 



m 



K* 



{tub + nip) 
X {K6K>i - ytbV;, [04 + aio] } . 



(132) 



(133) 



.G 



M{B'^K*'p') = _ili^/^,„^{(e^.e;,)(ms + m^.)^r''*K) 



X 



-{^p-Pb){^k - Pb) 



>K* 



(m^ + niK* 



i^lJ.iyaP^p^'KPBPK 



2V''-''' {ml) 
{niB + niK*) 



— i- 



VubV*sa2 - VtiMt*, X -(ag + 07) 



fK*mK- {(cp • eK)(niB + nip)Af-''{mj^*) 
-{ep . pB){eK ■ p^)^J^Z>M. _ ,V ..M 



(ms + mp) ^''"^ K p B p _^ 



04 — 2^10 



(134) 



M{B- 



K* 



G 

-i^fp'^p{(^o ■ e-){niB + mj^.)Af^^*(mJ) 



2A^- 



f a P 
l^p.uajS^O^-PBP- 



X 



3 ^ 

K&K>2 - VtbVZ X -(ag + 07) I 

i^fK*niK* {(eo • e_)(mB + mp)Af-''{mj^*) 
-(eo •Ps)(e- •Ps)-73: — ^"17^^ ^W/3^-^oPbPo 



21^ 



(m2) 



(ms + mi<:«) 



(135) 



(ms + nip) 
X {K&K>i - VtbVZ [04 + aio] } . 



{niB + nip) 
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X 



-(eo • PB)(e^ ■ Pb)—, ■ ^ «W/3eoe>sPo 



m 



K* 



X Wtt[-«4 + ^Oio] 



G 



{K6K>2 - Vt,V:,[2{a^ + as) + ^(09 + 07)} . (136) 



M{B-^K*-u) = -i^/K*m;^*{(eo-e_)(mB + m,)Af^'^(m|0 

[niB + m^) [mB + rriu,) 

X { VubV:,ai - VtbV* [a4 + aio] } 

2 



,2Af^^*(m2) . u . „ /3 2y^-^-^*(m2)l 
- eo • e_ • Pb) , - ^eM-/3e^e-PSp^7 — 

X [VubV:,a2-VtbVZ[2{as + a,) + ^{ag + a7)Y (137) 

(3) Pure penguin processes: 

M{B-^p-K*') = z^/K*mK*{(e,-e;,)(mB + m,)Af^''(m^,) 



i^p ■ PB){eK ■ Pb) 



(niB + nip) 

^W/34^»g . J^:^ [ V^t^V*: {«4 - -aio I . (138) 



(iriB + rup) 



, w ,2Af-(m|) 
-(e,-PB)(e.-pB) + 

. . g , 2V-^-(m|) | 
'"'^'^"'^'^''^^^^'^ {tub + m.) j 

X ^tbVjdjaa + as- ^(07 + 09) }• (139) 
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-{e,-PB){e,-PB) ^^^^^^^ 

..... r I, 

X 



jag + as - ^(ar + a9)| . (140) 
M{B- ^ p-(j)) = -\/2M(B° ^pV)- (141) 



M(5~ ^ ir*-0) = 

-(e,-p.)(6^-p.) 

-ie,..pe^e^KPBPK (^^ + ^^,)| 
X '^tft^ts |«3 + 04 + 05 - ^(07 + Og + aio)| . (142) 



M(S- ^ X*-X*°) = 

-ie„uai3eie2PBP2 —/ , r f ^tbVtd \ - -aio \ ■ (143) 
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